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ABSTRACT 

In this paper we present computations of the integrated metal-ion column densities produced in the 
post-shock cooling layers behind fast, radiative shock-waves. For this purpose, we have constructed 
a new shock code that calculates the non-equilibrium ionization and cooling; follows the radiative 
transfer of the shock self-radiation through the post-shock cooling layers; takes into account the re- 
sulting photoionization and heating rates; follows the dynamics of the cooling gas; and self-consistently 
computes the initial photoionization state of the precursor gas. We discuss the shock structure and 
emitted radiation, and study the dependence on the shock velocity, magnetic field, and gas metallicity. 
We present a complete set of integrated post-shock and precursor metal-ion column densities of all 
ionization stages of the elements H, He, C, N, O, Ne, Mg, Si, S, and Fe, for shocks with velocities of 
600 and ~ 2000 km s~^, corresponding to initial post-shock temperatures of 5 x 10^ and 5 x 10^ K, 
cooling down to 1000 K. We consider shocks in which the magnetic field is negligible [B = 0) so 
that the cooling occurs at approximately constant pressure ( "isobaric" ) , and shocks in which the 
magnetic pressure dominates everywhere such that the cooling occurs at constant density (isochoric). 
We present results for gas metallicities Z ranging from 10""^ to twice the solar abundance of heavy 
elements, and we study how the observational signatures of fast radiative shocks depend on Z. We 
present our numerical results in convenient online figures and tables. 

Subject headings: ISM:general - atomic processes - plasmas - absorption lines - intergalactic medium 
- shock waves 



1. INTRODUCTION 

Shock waves are a common phenomenon in astro- 
physics, and have a profound impact on the energetics 
and distribution of gas in the interstellar and intergalac- 
tic medium (McKec & HoUenbach 1980; Draine & Mc- 
Kee 1993 and references therein). A detailed understand- 
ing of the physical properties and observational signa- 
tures of shock-heated gas is needed, for example, in the 
study of diffuse gas in supernova remnants (e.g. ShuU 
& McKee 1979; Dopita et al. 1984); young stellar ob- 
jects (e.g. Dopita 1978); high velocity clouds (e.g. Fox et 
al. 2005); active galactic nuclei (e.g. Dopita & Suther- 
land 1995) - mostly in the context of narrow emission 
line regions, low ionization emission line regions, and ra- 
dio galaxies; the shock-heated "warm/hot" intergalactic 
gas (e.g. Savage et al. 2005; Furlanetto et al. 2005; Ccn 
& Ostriker 2006; Bertone et al. 2008); and in the virial 
and accretion shocks that may form in hot and cold-flow 
accretion into dark-matter halos and proto-disks in high 
redshift galaxies (e.g. Birnboim & Dekel 2003; Keres et 
al. 2005; Dekel et al. 2008). 

In this paper we present new computations of the non- 
equilibrium metal ionization states and integrated metal 
column densities that are produced behind fast, steady, 
radiative shock waves, and of the equilibrium columns 
produced in their radiative precursors. Our work is mo- 
tivated by recent Hubble Space Telescope (HST), Far Ul- 
traviolet Spectroscopic Explorer (FUSE), Chandra, and 

Electronic address: orlyg@tapir.caltech.edu 

^ School of Physics and Astronomy and the Wise Observatory, 
The Beverly and Raymond Sackler Faculty of Exact Sciences, Tel 
Aviv University, Tel Aviv 69978, Israel 

^ Theoretical Astrophysics, California Institute of Technology, 
MC 130-33, Pasadena, CA 91125 

^ Chandra Fellow 



XMM-Newton, detections and observations of hot gas 
around the Galaxy (CoUins et al. 2005; Fang et al. 2006), 
and of gas that may be a part of a 10^-10^ K "warm/hot 
intergalactic medium" (WHIM; Tripp et al. 2000, 2006, 
2007; ShuU et al. 2003; Richter et al. 2004; Sembach et 
al. 2004; Soltan et al. 2005; Nicastro et al. 2005; Savage 
et al. 2005; Rasmussen et al. 2006; Wilhams et al. 2006). 
An intergalactic shock heated WHIM is expected to be 
a major reservoir of baryons in the low-redshift universe 
(Cen & Ostriker 1999, 2006; Dave et al. 2001; Furlan- 
etto & Loeb 2004; Bertone et al. 2008; Danforth & ShuU 
2008). Important ions for the detection of such gas in- 
clude C IV, N V, O VI, O VII, O VIII, Si III, Si IV, 
Ne VIII, and Ne IX. 

Recent observations have confirmed the existence of 
warm/hot gas, both in the local Universe, and in more 
distant environments (e.g. Tripp et al. 2007; Savage et 
al. 2005). For example. Savage et al. (2005) reported 
on an absorption-line system at a redshift of z = 0.207, 
containing C III, O III, O IV, O VI, N III, Ne VIII, 
Si III, and Si VI. They invoked an equilibrium two- 
phase model, in which most of the ions are created in 
a warm (~ 2 x 10'* K) photoionized cloud, while O VI 
and Ne VIII are produced in a hotter (~ 5 x 10^ K), 
coUisionally ionized phase. The inferred temperature 
for the coUisional phase is within the range of tempera- 
tures where departures from equilibrium are expected to 
be important (Gnat & Sternberg 2007; hereafter GS07). 
While observations have confirmed the existence of the 
warm/hot clouds, detailed theoretical modeling of the 
physical properties in the gas is needed to determine 
what fraction of baryons they harbor, and whether they 
confirm the theoretical predictions regarding the shock- 
heated WHIM (Furlanetto et al. 2005). Investigating 
how the absorption line properties of fast radiative shock 
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waves depend on gas metallicity is particularly interest- 
ing in the context of intergalactic shocks, as some of the 
WHIM absorbers are expected to have significantly sub- 
solar metallicities (e.g. Gen & Ostriker 2006). 

We focus on fast shocks, with velocities, Vs , of 600 and 
2000 km s~^, corresponding to post-shock temperatures, 
Ts, of 5 X 10^ to 5 X 10^ K. After being heated by the 
shock, the gas emits energetic radiation which may later 
be absorbed by cooler gas further downstream, and by 
the unperturbed gas that is approaching the shock front. 
The shock self-radiation significantly affects the ioniza- 
tion states and thermal properties of the gas. 

In following the time-dependent ion fractions, we rely 
on the code and work presented in GS07 where we fo- 
cused on pure radiative cooling and ionization of gas 
clouds, with no external source of heating or photoion- 
ization. Below ^ 10^ K, cooling can become rapid com- 
pared to ion-electron recombination. When recombina- 
tion lags behind cooling the gas at any temperature tends 
to remain over-ionized compared to gas in collisional 
ionization equilibrium (Kafatos et al. 1973; Shapiro & 
Moore 1976; Edgar & Chevalier 1986; Schmutzler & 
Tscharnuter 1993; Sutherland & Dopita' 1993). In GS07 
we presented results for metal ionization states assuming 
collisional ionization equilibrium, as well as for the non- 
equilibrium ionization states that occur in radiatively 
cooling gas in which recombination lags behind cool- 
ing. Here we extend our work to compute the integrated 
metal-ion column-densities produced in steady flows of 
gas, in the post-shock cooling layers behind fast shocks. 
The computations we present in this paper build upon 
the results we presented in GS07. In addition to calculat- 
ing the non-equilibrium ionization and cooling, we follow 
the dynamics of the cooling gas and the radiative transfer 
of the shock self-radiation through the post-shock cool- 
ing layers. We include the effects of photoionization and 
heating by the shock radiation, and also self-consistently 
compute the initial photoionization state of the precursor 
gas. 

The diagnostic characteristics of gas cooling behind 
steady shock waves have been studied by many au- 
thors dating back to Cox (1972), Dopita (1976; 1977) 
and Raymond (1979). These works computed the time- 
dependent ionization and cooling in the post-shock cool- 
ing layers. ShuU & McKee (1979) added a self-consistent 
computation of the ionization states in the radiative 
precursor produced by shock radiation emitted in the 
up-stream direction. All of these investigations con- 
sidered intermediate or low-velocity shocks, with Vs < 
200 km s~^. Kang & Shapiro (1992) considered shocks 
with Vs up to 400 km s~^, for gas with primordial com- 
positions. For slow shocks, the neutral fraction in the gas 
entering the shock may be considerable, and this signifi- 
cantly affects the thermal and dynamical evolution of the 
post-shock gas. The metal-line signatures of slow shocks 
differ significantly from those of faster shocks, both be- 
cause collisional ionization and excitation behind cooler 
shocks are less efficient, and because the photoionizing 
self-radiation produced by the shocked gas is much less 
intense and less energetic. 

The observational signatures of fast shock waves 
{vs > 400 km s~^) were initially studied by Daltabuit, 
MacAlpine & Gox (1978), and Binette et al. (1985) as- 
suming collisional ionization equilibrium, later hy Dopita 



& Sutherland (1995; 1996), and more recently by Allen 
et al. (2008), who explicitly followed the time-dependent 
ion fractions, and took into account the photoionization 
and heating of both the "downstream" gas and radia- 
tive precursor by the shock self-radiation. Allen et al. 
present results for shocks with velocities between 100 and 
1000 km s^^, and with magnetic field parameters Bj^n 
between 10~^ and 10 /xG cm^/^. They considered five dif- 
ferent heavy element abundance-sets, all with a roughly 
solar metals-to-hydrogen ratios. They presented results 
for the integrated metal column densities and resulting 
emission line properties. 

None of the previous discussions investigated how the 
cooling column densities depend on the overall metallic- 
ity of the shocked gas. For the 600 and ^ 2000 km s~^ 
shocks that we consider, we present a complete set of re- 
sults for gas metallicities Z ranging from lO^"* to twice 
the solar abundance of heavy elements. We investigate 
how the observational signatures of fast radiative shocks 
depend on Z . Our emphasis is on the absorbing metal 
column densities that are produced in the post-shock 
cooling layers. We present computations for shocks in 
which the magnetic field is negligible [B = 0) so that 
the cooling occurs at approximately constant pressure 
(isobaric), and for shocks in which the magnetic field 
dominates the pressure everywhere so that the gas cools 
at constant density (isochoric). There are two key dif- 
ferences between isobaric and isochoric shocks that we 
discuss in detail in this paper. First, in isobaric flows 
the cooling times are shortened in the downstream gas 
due to the higher densities that obtain as the gas cools 
and compresses. This leads to a large suppression of the 
integrated columns densities of species that are produced 
mainly in the cooler parts of the flows, in isobaric versus 
isochoric shocks. Second, the effects of photoionization 
by the shock radiation are much more important in iso- 
choric shocks where the "ionization parameters" remain 
large because the gas densities do not grow as the ra- 
diation is absorbed. Photoionization in isochoric shocks 
leads to increased ionization in the cooler downstream 
absorbing layers, over and above what is produced by 
recombination lags in the rapidly cooling gas. The mag- 
nitude of these effects is sensitive to the metallicity, via 
its control of the cooling rates and the associated shock 
self-radiation. 

In our computations we make several simplifying as- 
sumptions. First, we assume steady state, one dimen- 
sional, plane-parallel shocks. Fast radiative shocks are 
known to be subject to thermal instabilities, both lo- 
cally - where small inhomogeneities in the gas grow to 
produce clumps and create secondary shocks, and glob- 
ally where the shock velocity oscillates, thus affect- 
ing the dynamics of the flowing gas (e.g. Draine & Mc- 
Kee 1993). Global shock instabilities are expected to oc- 
cur for shock velocities in excess of 140 km s~^, based on 
both analytical (e.g. Chevalier & Imamura 1982; Innes et 
al. 1987; Innes 1992; Toth & Draine 1993) and numerical 
work (e.g. Langer et al. 1981; Gaetz et al. 1988; Toth & 
Draine 1993; Strickland & Blondin 1995; Sutherland et 
al. 2003; Pittard et al. 2005). Global oscillations of the 
shock velocity may affect the absorption and emission 
line signatures of the post-shock cooling gas, and it is 
unclear to what extent predictions based of steady-state 
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shock models represent the complexity of real physical 
shocks. However, as noted by Allen et al. (2008), full 
three-dimensional models of unsteady and dynamically 
evolving shocks, including radiative transfer and time- 
dependent evolution of the ion fractions in the post-shock 
cooling layers is computationally very complex. This 
is beyond the scope of this work. Here we construct 
steady-state one-dimensional models as approximations. 
In our steady state, one-dimensional models, we self- 
consistently take into account the time dependent ion- 
ization, cooling, radiative transfer, and dynamics of the 
post-shock gas. 

Second, we assume that the electron and ion tem- 
peratures are equal throughout the flow. Recent ob- 
servational evidence has suggested that in "Balmer- 
dominated" shocks with ti^ > 400 km the electrons 
and ions emerge from the shock with different post- 
shock temperatures (e.g. Ghavamian et al. 2007). In 
600 and 2000 km s~^ shocks, Ghavamian et al. esti- 
mate that the post-shock electron temperatures are ~ 0.5 
and 0.05 the proton post-shock temperatures. The va- 
lidity of the standard electron-ion equipartition assump- 
tion has therefore been questioned (Ghevalier et al. 1980; 
Ghavamian et al. 2001, 2007; Yoshida et al. 2005; Heng 
& McCray 2007; Heng et al. 2007; Rakowski et al. 2008; 
Raymond et al. 2008). However, the electron-proton 

equipartition time t^p — 8.4x lO'^Tg/^/rie years (where Ue 
is the electron density, and T^r — Te/IQ^ K, Spitzer 1962; 
Yoshida et al. 2005), is much shorter than the cooling 

times, tc w 3.9 x lO^Tj^'^ /{rieZ) (GS07), for the physical 
conditions that we consider, and any initial temperature 
difference will be rapidly removed. The electron and ion 
velocity distributions are therefore expected to remain 
closely coupled during the cooling process. For example, 
for solar metallicity {Z = 1), we estimate, conservatively, 
that for 600 and 2000 km s~^ shocks the initial ratios of 
the equipartition to cooling times are 4.3 x 10"'^ and 
7.2 X 10~^ respectively"*. The ratio tgp/tc is smaller for 
lower Z. We therefore assume Tg = 21on everywhere. 

Third, we assume that any initial dust is rapidly 
destroyed (e.g. by thermal sputtering) after passing 
through the shock front, on a time scale that is much 
shorter than the cooling time. If a high dust mass can 
be maintained, cooling by gas-grain collisions may dom- 
inate the total cooling at high temperatures (Ostriker 
& Silk 1973; Draine 1981). Draine & McKee (1993) es- 
timate that as much as 15% of the initial thermal en- 
ergy may be radiated away before the grains are com- 
pletely destroyed. Even in this case, the integrated cool- 
ing column densities are only affected at a level of < 10% 
(except for the highest ionization state of each element 
which may be affected by < 30%). More recently, the 
work of Smith et al. (1996) has suggested that the dust 
sputtering destruction time-scale is shorter than the cool- 

In estimating t^p we set T = Ts- This gives an upper bound 
on the equipartition time. In estimating tc we set T equal to the 
initial electron temperature Te. This gives a lower bound on the 
cooling time. Thus, for example, for a 600 km shock, with 
= 5 X 10^ K, and for which initially Te = 0.5Tp (where Tp is 
the proton temperature), and with Ts = (Te -|- Tp)/2, it follows 
that Te = (2/3)Ts = 3.3 x 10"^ K. For these values of Ts and Te 
our conservative estimates for the equipartition and cooling times 
are tep — 3.0 X lO^/ne and tc = 7.1 X IQ^ /{ueZ) years. This gives 
tep/tc = 4.3 X 10"^ for Z = 1. 



ing time when T > 3 x 10^ K. We consider higher shock 
temperatures, and assume that any dust is instanta- 
neously destroyed. The post-shock gas then evolves with 
constant gas-phase elemental abundances as specified by 
the metallicity, Z . 

The outline of our paper is as follows. In §2 we write 
down the equations of ionization, dynamics, cooling and 
heating, and radiative transfer that we solve in our com- 
putations. We also discuss our treatment of the initial 
ionization states in the gas approaching the shock front, 
and describe our numerical method. In §3 we discuss the 
shock structure and emitted radiation, and investigate 
how they depend on the controlling parameters, includ- 
ing the gas metallicity, shock velocity, magnetic field, and 
gas density. In §4 we present our computations of the 
non-equilibrium ionization states in the post-shock cool- 
ing layers. We describe how photoionization by the shock 
self-radiation and departures from ionization equilibrium 
affect the ion fractions in the cooling gas. In §5 we discuss 
the radiative cooling and heating efficiencies. Finally, in 
§6 we present the integrated metal-ion column-densities 
that are produced in steady flows of cooling gas behind 
fast radiative shocks, for gas metallicities between 10~^ 
and 2 times solar abundances. We discuss how the col- 
umn densities depend on the gas metallicity, shock ve- 
locity, and magnetic field. In § 7 we compute the equi- 
librium metal-ion column densities that are produced in 
the photoionized precursors of these fast radiative shocks. 
We summarize in §8. 

In this work we have generated a large num- 
ber of tables containing our numerical results. 
These are a vailable as online data files at 
http://wise-obs.tau.ac.il/~orlyg/shocks/ 

2. BASIC EQUATIONS AND PROCESSES 

We are interested in studying the evolving ionization 
states in steady flows of cooling gas. We focus on the ion- 
ization states in post-shock cooling layers behind steady 
fast radiative shock waves. The gas is heated to some 
high temperature Ts > 5 x 10^ K by the shock, and 
then cools and recombines as it flows away from the 
shock front. In our models we follow a gas element 
as it advances through the post-shock flow. If the gas 
cools faster than it recombines, non-equilibrium effects 
become significant, and the gas remains over-ionized. 
Non-equilibrium ionization leads to a suppression of the 
cooling rates. As in GS07, we compute this coupled 
time-dependent evolution. We describe the ionization 
and cooling in § 12.11 and § 12.31 below. 

The dynamics of a gas parcel that moves through the 
flow is determined by the continuity equations for the 
mass-, momentum-, and energy-fluxes (see § 12. 2p . We 
explicitly follow these equations to derive the local gas 
velocity, density, temperature, and pressure. 

As gas enters the shock, it is heated to a high tem- 
perature Tg. The gas then flows away from the shock 
front, and gradually radiates its thermal energy. For fast 
shock waves, this radiation is intense and energetic. The 
shock self-radiation is later absorbed by cooler gas fur- 
ther downstream, providing a source of heating and pho- 
toionization, which greatly affects the shock structure 
and the ion fractions in the cooling gas. To compute 
the local ionization and heating rates everywhere in the 
flow, we follow the transfer of the shock self-radiation as 
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it advances downstream (see § I2.4|) . 

The radiation field that is emitted by the shock pene- 
trates the gas approaching the shock front, and creates a 
radiative precursor that is photoionized by this radiation 
(see § I2.5|) . The ionization states in the radiative pre- 
cursor are in photoionization equihbrium with the shock 
self-radiation, and set the ion fractions in the gas enter- 
ing the shock front. As we describe in § I2.6[ we iterate 
to find a self-consistent solution for the initial ionization 
state and the shock self-radiation. 

As we describe in § [TJ we assume that the shocks 
are steady, and that after passing the shock fronts the 
electron and ion temperatures are equal. In computing 
the cooling rates and abundances, we neglect cooling by 
molecules and dust. 

2.1. Ionization 

As in GS07, we consider all ionization stages of the ele- 
ments H, He, C, N, O, Ne, Mg, Si, S, and Fe. We include 
collisional ionization by thermal electrons, radiative re- 
combination, dielectronic recombination, and neutraliza- 
tion and ionization by charge transfer reactions with hy- 
drogen and helium atoms and ions (GS07, and references 
therein). Here we add statistical charge transfer rate 
coefficients for high-ions with charges greater than -1-4 
(Kingdon & Ferland 1996). In this work we also include 
photoionization (Verner et al. 1996) and multi-electron 
Auger ionization processes (Kaastra & Mewe 1993) in- 
duced by the X-ray photons emitted by the cooling gas. 
(As we discuss in §4.6, the effects of Auger processes are 
small.) 

The time-dependent equations for the ion abundance 
fractions, Xi, of element m in ionization stage i, at every 
point in the flow, are 

-Xi [{q.i + ai)nc + + kf^n-^o + kf^n^^o 
+kf^n-ti+ + kf^n-tic+] ■ 

In this expression, qi and are the temperature- 
dependent rate coefficients for collisional ionization and 
recombination (radiative plus dielectronic), and k^^, 

and k^^ are the rate coefficients for charge transfer 
reactions with hydrogen and helium that lead to ioniza- 
tion or neutralization. The quantities nno, nH+, n-n^o, 
riHc+i and Uc are the particle densities (cm~'^) for neu- 
tral hydrogen, ionized hydrogen, neutral helium, singly 
ionized helium, and electrons, respectively. Tj^i are the 
local rates of photoionization of ions j which result in the 
ejection of i—j electrons. Fj are the total photoionization 
rates of ions i due to externally incident radiation. 

For each element m, the ion fractions Xi = 
n-i^m/inHAjn) must at all times satisfy 

E^' = l ' (2) 

where n^^m is the density of ions i of element m, nn — 
njjo 4- fiH+ J and Am is the abundance of element m rel- 
ative to hydrogen. 

2.2. Dynamics 



As we discussed in § [1] we assume one-dimensional, 
steady-state radiative shocks. We neglect any instabil- 
ities that may form, even though fast radiative shocks 
are knows to be subject to a global oscillatory instability 
for shock velocities above ^ 140 km s~^ (e.g. Ghavamian 
et al. 2007). We also assume that the electron and ion 
temperatures are equal everywhere. As discussed in §1, 
the electron-proton equipartition times are much shorter 
than the cooling times for the shock velocities we con- 
sider. The electron- and ion-velocity distributions are 
therefore expected to remain closely coupled during the 
cooling process. 

Under these assumptions, we follow the standard 
Rankinc-Hugoniot conditions (e.g. Cox 1972): 

pv = const 

B/ p = const 

pv'^ + P + ^B'^ = const 

ipw"^ H — hjP''^ + Pv + j^B'^v + J^{nenHA — nT)dx — const 

^ . . .(3) 

where p is the mass density, v the velocity in the direc- 
tion of the flow, P is the gas pressure, B is the magnetic 
flcld perpendicular to the direction of the flow, n is the 
particle density, A is the cooling efficiency (erg s~^ cm'^), 
and T the heating rate (erg s^^). These equations rep- 
resent the conservation of the mass-, momentum-, and 
energy-fiux in the flow, and the condition of "frozen in" 
magnetic field. 

The assumption of a strong shock implies the famil- 
iar jump conditions (e.g. Draine & McKee 1993) relat- 
ing the pre-shock and post-shock physical properties: 
no — 'iripre, va = 1/4 Vs- In these expressions no and 
vq are the post-shock particle density and velocity, Upre 
is the pre-shock density, and Vg is the shock velocity in 
the frame of the pre-shocked gas. The shock temperature 
and velocity are related by Ts — Sp^v^ /16kB (e.g. McKee 
& HoUenbach 1980), where p is the mean mass per par- 
ticle, and fee is the Boltzmann constant. 

We consider two limiting cases. First, we consider flows 
in which B = everywhere. The gas then cools at ap- 
proximately constant pressure (see § 12. 3p . Second, we 
study the case where the magnetic field is "dynamically 
dominant" everywhere, with B j ^fp ^ tig, so that it dom- 
inates the pressure throughout the flow (see equation [3]). 
In this case, the condition of "frozen-in" magnetic fleld 
implies constant density evolution. 

Realistic shocks may have some intermediate value for 
the magnetic fleld, and the observed signatures of such 
shocks will lie between the two limiting cases presented 
in this paper (e.g. Dopita & Sutherland 1996; Allen et 
al. 2008). 

2.3. Cooling 

The ionization equations are coupled to an energy 
equation for the time-dependent heating and cooling, and 
resulting temperature variation. 

We follow the electron cooling efficiency, A(T, a;^,^) 
(erg s~^ cm'^), which depends on the gas temperature, 
the ionization state, and the total abundances of the 
heavy elements specified by the metallicity Z. As in 
GS07, we adopt the elemental abundances reported by 
Asplund et al. (2005) for the photosphere of the Sun, 
and the enhanced Ne abundance recommended by Drake 
& Testa (2005). We list our assumed solar abundances 
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TABLE 1 
Solar Elemental 
Abundances 



Element 


Abundance 




(X/H)0 




Carbon 


2.45 X 10" 


'4 


Nitrogen 


6.03 X 10" 


-5 


Oxygen 


4.57 X 10" 


-4 


Neon 


1.95 X 10" 


-4 


Magnesium 


3.39 X 10- 


-5 


Silicon 


3.24 X 10" 


■5 


Sulfur 


1.38 X 10" 


■5 


Iron 


2.82 X 10" 


-5 



in Table [T] In all computations we assume a primordial 
helium abundance Ane = 1/12 (Ballantyne et al. 2000), 
independent of Z . 

The electron cooling efficiency includes the removal 
of electron kinetic energy via recombinations with ions, 
collisional ionizations, coUisional excitations followed by 
prompt line emissions, and thermal bremsstrahlung. We 
do not include the ionization potential energies as part of 
the internal energy but instead follow the loss and gain 
of the electron kinetic energy only (see GS07). 

We also follow the heating rate, T{xi,Z, J^) (erg s~^) 
due to absorption of the shock self-radiation, J^, by gas 
further downstream. The net local cooling rate per vol- 
ume is given by n^nH^ — nT , where n is the total gas 
density. 

For an ideal gas, the pressure P — nk^T, and the 
thermal energy density u = 3/2nfcBT, where T is the 
gas temperature, and n is the total particle density. We 
follow the gas pressure, density, temperature and velocity 
as is flows away from the shock front. 

When the magnetic field is set to zero, the Rankine- 
Hugoniot conditions relate the net local cooling to the 
gas deceleration: 



3c^ 
2 



V dv 



dt 



P 



(4) 



(e.g. Shu 1992), where = (5/3) P/p is the sound 
speed. We use equation ^ to follow the evolution of the 
velocity along the flow. We later use the mass and mo- 
mentum continuity conditions to derive the gas density 
and pressure via. 



p = pqvq/v and 
P = /3owo(4wo - v), 



(5) 



where p — np and v are the local gas mass density and 
velocity, po and vq are the post-shock mass-density and 
velocity, and p is the mean mass per particle. Right after 
passing the shock front, the gas density is determined by 
the jump conditions (po — ^Ppre] vq = Vs/'^)- The gas 
pressure is then given by Pq = 'ipov^ (see Equation [5]). 
As the gas radiates all of its internal energy, its veloc- 
ity becomes very small, and the final pressure is then 
Poo = 4poWo — 4/3 Po- Thus, B = Q corresponds to 
approximately isobaric dynamics. 

In the limit of a dynamically dominant magnetic field 
(i?/y/p S> Vs), the Rankine-Hugoniot conditions imply 
constant density and velocity throughout the flow, so 
that p — pq and v = vq everywhere. This is an isochoric 



flow. We refer to this as the "strong-B" limit. Energy 
conservation then implies that the radiated energy equals 
the loss of thermal energy, and relates the net local cool- 
ing to the pressure change by 



3dP 
2~dt 



(6) 



(see equation [7] in GS07). 

Given a set of nonequilibrium ion abundances, Xi{T), 
and the gas metallicity Z, we use the cooling and heating 
functions included in Cloudy (ver. 07.02.00; Ferland et 
al. 1998) to calculate A{T,Xi,Z) and T{xi,Z,J^). 

2.4. Radiative Transfer 
We follow the radiative transfer equation, 

dl„{p) 



dr,. 



Iv{p) - Sp, 



(7) 



intensity (erg s ^ cm ^ Hz 
source function, 



to evaluate the intensity of the radiation along the flow. 
In equation ([7]), Iv{p) is the direction-dependent speciflc 

sr"-'^). Si, = (.yjoLy is the 
is the specific emissivity (or emis- 
sion coefficient, erg s"^ cm"-^ Hz^^ sr^^), is the ab- 
sorption coefficient (cm~^), Ty is the optical depth, and 
p = cos(ff). 

We use a Gaussian quadrature scheme to evaluate the 
local intensity of the radiation as a function of distance 
from the shock front (Chandrasekhar 1960). We follow 
the specific intensity l^i^p) along ten downstream direc- 
tions between p= \ (parallel to the shock velocity) and 
/i = (perpendicular to the velocity). We divide the 
flow into thin slabs of thickness Z(T) (see § 12. 6p . For 
each slab, we use Cloudy to evaluate the absorption (ay) 
and emission (ty) appropriate for the local conditions 
that we calculate. We use our local non-equilibrium ion- 
ization states and temperature at the upstream edge of 
each slab as input for the Cloudy models. The Cloudy 
output includes the continuum specific emissivities, a list 
of line emissivities, and the continuum absorption coeffi- 
cients as a function of energy^. 

Once the opacities and emissivities for a slab are 
known, we compute the radiative transfer for each of 
the 10 values of p. The optical depth in the downstream 
direction is given by r° = / aj^, where / is the thick- 
ness of the slab in the downstream direction. The opti- 
cal depth along the different directions is then given by 
Ty{p) = tII p. The source function Sy = tyjay. Given 
an input specific intensity, /™(/z), the specific intensity 
at the downstream edge of a slab is 



\t^) 



ili:\p) ^ Sy) + Sy 



(8) 



/°"* is then used as input radiation for the next slab. 
Finally, the mean intensity, Jy = ^ Iy{p)dVl, is com- 
puted locally, and used in evaluating the local heating 
and photoionization rates. 

2.5. Radiative Precursor 

The initial conditions for the shock depend on the ion- 
ization states in the radiative precursor. The precursor 

^ It does not include line opacities. In the Cloudy computations, 
we therefore do not include continuum pumping of the lines, as 
these would not be consistently absorbed. 
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is photoionized by the intense radiation field created by 
the shocked gas (Shull & McKee 1979). If the velocity of 
the ionization front in the gas approaching the shock, 

1 f°° Ju 

(where Upre is the preshock density) is larger than the 
shock velocity Vs, a stable photoionization equilibrium 
radiative precursor will form. This condition is met for 
shock velocities > 175 km s~^ (Dopita & Sutherland 
1996), and a stable radiative precursor will therefore form 
for the shock velocities that we consider here. 

We self consistently calculate the ionization states in 
the radiative precursor (see § [7]) . Since the photoionizing 
radiation is emitted by the shocked gas, iterations are 
required to obtain a self consistent solution. In the first 
iteration, we assume that the gas starts out in collisional 
ionization equilibrium (CIE) at the shock temperature, 
Tg. We then compute the resulting shock model, and 
follow the radiative transfer in the downstream direction 
(see §1231) • 

The self-radiation builds up with distance from the 
shock front, as more and more emitting gas contributes to 
the intensity. However, at some distance from the shock 
front, the gas becomes optically thick, and the intensity 
begins to decrease. We label the distance at which the 
gas becomes thick at the Lyman limit hhick- At hhick, 
most of the initial energy has been emitted, but it has not 
yet been absorbed. We assume that the mean intensity 
at this point, Ju{hhick), is similar to the intensity enter- 
ing the radiative precursor, and use Jv(hhick) to compute 
the ionization states in the radiative precursor. 

Since our shock velocities are high enough to ensure the 
formation of a steady equilibrium photoionized precur- 
sor, we use the photoionization code Cloudy to compute 
the ion fractions in the radiative precursor. 

We then use these ionization states as initial condi- 
tions for a second iteration. This process can be repeated 
until the resulting photoionizing radiation Ju{hhick) has 
converged. We find that the second iteration is suffi- 
cient, and that although the ionization states in the pre- 
cursor are significantly underionized relative to CIE at 
the shock temperature, the radiation field created by the 
shocked gas is not substantially altered. The shock struc- 
ture (e.g. the temperature as a function of time, and the 
integrated column densities) is also similar in the first 
and second iterations. 

The assumption that the upstream radiation field 
equals Ju{khick) is based on the fact that the flowing 
gas is optically thin between the shock front and hhick- 
Nevertheless, geometrical effects as well as scatterings in 
the downstream gas, may increase the intensity of the 
radiation in upstream directions. Our assumption may 
only lead to an underestimate of the intensity of radiation 
photoionizing the precursor. The true ionization states 
in the precursor may therefore be higher than we com- 
pute, but still lower than the CIE values used in the first 
iteration. Since the two iterations lead to similar shock 
structures and mean intensities, we conclude that this 
assumption does not significantly affect the solution. It 
may have a minor impact on the ionization states imme- 
diately following the shock-front before the gas adjusts to 
CIE at the shock temperature. This will have a negligible 
impact on the integrated column densities (see § ^ . 



2.6. Numerical Method 

The abundance equations ([T]) and the fiow equation 
(eq. [4] when _B = or eq. [6] in the strong-i? limit) 
are a set of 103 coupled ordinary differential equations. 
When _B = 0, we advance the numerical solution in small 
velocity steps Av = ev, where e < 0.005, and v is the 
velocity associated with a temperature T{v). For strong- 
B isochoric conditions we advance the solution in small 
pressure steps AP = eP, where P is the gas pressure 
associated with a temperature T{P). 

For any temperature T, we compute the total cooling 
and heating rates by passing the current non-equilibrium 
ion fractions Xi{T), temperature T, and mean intensity 
Jjy to the Cloudy cooling and heating functions. We 
then compute the time interval At associated with the 
velocity change Av = ev (when i? = 0) or pressure 
change AP = eP (for strong- _B). In some cases, ad- 
ditional constraints were set when determining the time 
step, to ensure numerical accuracy and computational ef- 
ficiency. As the gas approaches thermal equilibrium, the 
time steps At become very large, and may exceed the re- 
combination time. We therefore apply an upper limit on 
the time steps, Aimax = max(0.05 x t,10^-^ /uq/Z), that 
depends on the current time, and on the gas density and 
metallicity. We further demand that each time step will 
be at most 3 times the previous step. We find that this 
choice prevents significant numerical "noise" , but allows 
the computation to proceed efficiently. 

We integrate equations ([T]) over the interval At using a 
Livermore ODE solver (Hindmarsh 1983, see GS07). We 
assume that over the time step At the velocity (or pres- 
sure) evolves linearly with time. In the integration, the 
estimated local errors on the fractional ion abundances 
are controlled so as to be smaller than 10"^, 10^^, and 
10^^, for hydrogen, helium and metals, respectively. 

The radiation intensities in the different directions are 
stored as vectors indicating the intensities at different en- 
ergies. The energy grid contains 565 data point between 
an energy of 0.001 Ryd and ^ 850 Ryd. The energy grid 
is "dense" near ionization edges to ensure an accurate 
computation of the photo-absorption rates. 

When computing the ionization and heating rates, we 
assume that the intensity of the radiation field is constant 
within a slab of thickness I = v x At, and is equal to the 
intensity of the radiation entering that slab. For each 
time step, we use the ion fractions at the beginning of 
the slab to compute the emissivities and opacities of the 
gas in the slab. In following the radiative transfer, we 
assume that the emissivities and opacities are constant 
within the slab. 

In the first iteration, we assume that the gas starts at 
collisional ionization equilibrium at the shock tempera- 
ture. We then compute the self radiation emitted by the 
cooling gas. We use this radiation field to compute the 
initial photoionization equilibrium ion fractions for the 
second iteration (see § 12. 5|) . The equilibrium photoion- 
ization solution is significantly underionized relative to 
CIE at the shock temperature. Very rapid changes there- 
fore occur as the gas adjusts to CIE just below the shock 
temperature. In order to follow this rapid evolution ac- 
curately, we set e — 0.0001(-^^J2_)i-5 during this rapidly- 
evolving period. We find that the mean intensity emitted 
by the cooling gas converges by the second iteration to a 
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level of 10%, and no further iterations are required. 

3. SHOCK STRUCTURE AND SCALING 
RELATIONS 

In this section we describe the shock structure, and 
how it depends on the controlling parameters: shock 
temperature (or velocity), gas metallicity, and magnetic 
field. We show results for two values of shock temper- 
atures: 5 X 10^ K {vs ~ 600 km s^), and 5 x 10^ K 
{vs — 1920 km s""'^). We explore five different values of 
the gas metallicity Z, from 10"'^ to 2 times the metal 
abundance of the Sun. For each shock temperature and 
gas metallicity we study the shock structure and ion frac- 
tions for the B — 0, and strong-S limits. 

Figure [1] shows the temperature profiles of the post 
shock cooling layers for the different cases that we study. 
Panel (a) shows results for = 5 x 10^ K, for a "strong- 
s'' isochoric shock. The different curves show results 
for different gas metallicities. The horizontal axis shows 
no X t - the initial post-shock hydrogen density times 
time. This scheme makes the results nearly independent 
of density, as we discuss in section 3.5 below. 

3.1. A strong-B (isochoric), = 5 x 10^ K, shock at 
a metallicity of Z = 2. 

To illustrate the behavior of the post-shock cooling gas 
in the strong-_B limit, we first focus on the results for 
T, = 5 X 10^ K and Z = 2. 

The initial ionization states of the gas that enters the 
shock, are set by photoionization equilibrium of the pre- 
cursor gas with the shock self radiation. This gas has 
an equilibrium temperature which is significantly lower 
than the shock temperature. 

As the precursor gas passes through the shock front, its 
temperature abruptly rises to Tj, leaving the gas under- 
ionized relative to CIE at the shock temperature. The 
gas then very rapidly adjusts to CIE at a temperature 
close to Tg. During this phase, the gas radiates very effi- 
ciently, as the hot thermal electrons efficiently excite the 
low-energy transitions of the under-ionized gas. Dopita 
and Sutherland (1996) refer to this phase as the "ioniza- 
tion zone" . The cooling efficiencies are about 60 times 
higher than at CIE. The evolution to CIE is very rapid, 
and occurs within no x t = 10^^ cm~'^ s. After this time, 
the ion fractions and cooling efficiencies reach CIE, at a 
temperature very close to Tg (> O.OSTg). 

For the shock temperatures that we consider (> 5 x 
10'^ K), the CIE cooling efficiency at ~ is low. The 
gas therefore stays hot for a long time. This phase is 
the "hot radiative zone" (e.g. Draine & McKee 1993; 
Dopita & Sutherland 1996). Metal line emissions (and 
bremsstrahlung emission for low gas metallicity) domi- 
nate the cooling in this zone (Sutherland & Dopita 1993; 
GS07; see § ID). 

The hot radiative phase ends once the gas cools to a 
temperature for which the cooling efficiency is high, T < 
10^ K. The shocked gas therefore emits most of its initial 
thermal energy while in the hot radiative phase. Figure[2] 
shows the mean intensity of the radiation field emitted 
by the shocked gas. The gray curve shows results for 
Z ^ 2 times solar metallicity gas. For Z = 2, the relative 
contribution of resonance lines is very large, and the line 
to continuum contrast is high. The spectrum shows a 
very prominent "UV-bump" created by numerous UV 



emission lines, with excess radiation between ~ 1—5 Ryd. 
The line contribution is also significant at far-UV and 
even X-ray energies. This radiation is absorbed by cooler 
gas further downstream, providing a source of heating 
and photoionization. 

Once the temperature drops sufficiently to bring the 
gas closer to the cooling-efficiency peak, the tempera- 
ture dechne becomes very rapid, and the gas cools to 
a temperature of a few x 10'' K. If cooling becomes 
faster than recombination, departures from equilibrium 
occur, and the gas tends to stay over-ionized (Suther- 
land & Dopita 1993; GS07). We refer to this phase 
as the "non-equilibrium cooling zone" (e.g. Dopita & 
Sutherland 1996). During this rapid cooling stage, metal 
line emissions (and Hydrogen-Helium line emission for 
Z < 10^^) dominate the cooling. 

In Figure [3] we display various physical parameters in 
the cooling gas. The upper panel shows the tempera- 
ture as a function of (linear) time, starting at the non- 
equilibrium cooling zone. As the gas cools, hydrogen 
starts to recombine (see panel c). Eventually, the neu- 
tral hydrogen fraction is high enough that it allows for 
efficient absorption of the shock radiation. This occurs 
when xh i ~ lO"'^. This significantly raises the heating 
rate, as is shown by the gray line in panel (e). The in- 
tegrated H I column density measured from the shock 
front is shown in panel (d). 

Once the heating rate reaches the cooling rate (see 
panel e), the rapid temperature decrease stops, and the 
gas enters a "plateau" in which the temperature remains 
roughly constant, and the radiation is gradually ab- 
sorbed. This is the "photoabsorption plateau" (e.g. Do- 
pita & Sutherland 1996; also called the "recombination 
zone" [ShuU & McKee 1979] or the "thermahzation zone" 
[Draine & McKee 1993]). As the radiation is removed by 
absorption, the neutral fraction of the gas rises. In this 
optically thin and ionized part of the plateau, the heat- 
ing rate changes very slowly^. Once the neutral fraction 
is of order unity (> 0.1), the optical depth increases and 
the heating rate declines more rapidly. We refer to the 
part of the plateau in which hydrogen is still ionized as 
the warm ionized medium (WIM) plateau. 

The depth of the "WIM plateau" depends on the flux, 
i^jy, of ionizing radiation entering the plateau. The thick- 
ness of the ionized region will be 

^wiM = / diy, (10) 

HHTlea J hv 

where a is the recombination coefficient. The time spent 
in the WIM plateau is then ^vvim = ^wim/'^o — 2.5 x 
lO^"* s, following a 5 X 10^ K, = vjA = 150 km s-\ 
no = 0.1 cm~'^ isochoric shock. 

For the set of parameters considered here (T^ = 5 x 
10^ K, Z = 2 solar, strong-!?), after the hydrogen re- 
combines the gas rapidly cools down to a temperature 
T\ow — 1000 K, at which we terminate the computation. 
For other (lower) gas metallicities, an additional warm 
neutral medium (WNM) plateau may form, with a tem- 
perature between 8000 and 10, 000 K. The depth of the 
"WNM plateau" depends on the gas X-ray opacity, which 

^ The heating rate T ~ x^^^a < E >, where < E > is the 
mean photoelectron energy, and a the recombination coefficient. 
While the gas is optically thin and jj ^ 1, the heating rate is 
close to constant. 
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Fig. 1. — Temperature versus post-shock hydrogen density times time (cm ^ s), for gas metallicities ranging from 10 to 2 times solar. 
The upper panels show results in the strong-S limit (isochoric) and the lower panels show results for B = (approximately isobaric). The 
left hand panels show results for = 5 X 10^ K, and the right hand panels for Ts = 5 X 10^ K. The results presented in this figure were 
computed assuming a post-shock hygrogen density, njj = 0.1 cm~^, but are nearly independent of density (see § 3.5). 




E (Ryd) 

Fig. 2. — Emitted mean intensity (erg s~^ cm~^ Hz~^) versus 
energy (Ryd) for Ts = 5 X 10^ K in the strong-B (isochoric) limit. 
The gray curve shows the emitted spectrum for Z = 2, the dark 
dashed curve for Z = 10~^, and the dark solid curve for Z = 10^''. 



is set by the gas metallicity. The WNM plateau is more 
extended for lower metallicity gas (see Figure [1]) . 

The various components in post-shock cooling lay- 
ers (e.g. Draine & McKee 1993; Dopita & Sutherland 
2003) are illustrated schematically in Figured Here we 
have made the further distinction between warm ionized 
(WIM) and warm neutral (WNM) plateaus in the down- 



stream absorbing layers. 

Throughout most of the photoabsorption plateau and 
the final cooling following it, the radiation field is ab- 
sorbed slowly enough that the gas has time to approach 
photoionization and thermal equilibrium. The photoab- 
sorption time scale is given by tabs — ^na->v ' ^^'^ 
recombination time scale by tree = As long as 

tree << tabs, the gas stays at a state of near photoion- 
ization equilibrium with the local radiation field. 

Once the temperature drops below ~ 10** K, Lya emis- 
sion and other allowed transitions become very ineffi- 
cient, and the gas cools primarily via fine-structure (FS) 
line emissions, including O III 88.3/im, O III 51.8/im, 
Ne III 15.5^m, Si II 34.8/.tm, S III 33.5/zm, and 
C II 157.7/im. (We do not include molecules or dust 
in this computation). 

As mentioned above, for Z = 2 metal line cooling dom- 
inates throughout. Metal lines are important even for 
r ~ 10^ K where Lya cooling dominates in GIF. The 
Lya cooling efficiency is suppressed because the hydro- 
gen is still mostly ionized, while the existence of over- 
ionized metal species, such as Ne''^"'", 0^+, S^+, allows for 
efficient line cooling that is not available in GIE. Initially, 
the heating in the plateau is due to H I photoionization 
(forrtoxt < 5.3x10^^ cm^"^ s), then He I photoionization 
(for 5.3 X 10^^ < no X t < 7.4 X 10^^ cm^^ s), and then 
by photoionization of metals, mostly by neutral oxygen 
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Fig. 3. — Shock structure versus time for Ts = 5 X 10 K, for 
Z = 2 solar, and for a post-shock hydrogen density no = 0.1 cm""^, 
in the strong-B Hmit. (a) Temperature (K). (b) Electron fraction, 
(c) Neutral hydrogen fraction, (d) Integrated H I column den- 
sity (cm~-^) from the shock front, (e) Local cooling (dark) and 
heating (gray) rates per volume (erg s~^ cm~^). 

(for no xt > 7.4 x 10^^ cm~'^ s). These transitions occur 
as the photo-absorptions remove lower energy photons 
from the radiation field. For photon energies greater than 
1.8 Ryd, helium absorbs more efficiently than hydrogen, 
and at yet higher energies oxygen is more efficient than 
helium. 

Below ~ 10"' K, the cooling rate per volume generally 
decreases due to the decreasing electron density. Below 
^ 6000 K, the cooling rate per electron (n/f A) is constant 
to within 30% (for Z = 2), with small variations due 
to recombinations that increases the abundance of the 
dominant coolants (e.g. C II 157.7/^m and Si II 34.8 
mum). 

Because the radiation field controls both the heating 
rate, and the electron fraction (which influences the cool- 
ing rates), heating and cooling remain strongly coupled 
for the remainder of the cooling process. While both the 
heating and cooling rates decrease with time, they re- 
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Fig. 4. — A schematic plot of shock temperature versus time. 
The different zones are shown (from left to right): The preshocked 
upstream gas is initially in a photoionized "radiative precursor"; 
After passing the shock front, the gas enters the "Hot Radiative" 
zone, in which cooling is slow. It later passes through the "Non- 
Equilibrium Cooling" zone, in which cooling is rapid, and depar- 
tures from equilibrium ionization and cooling may occur. After 
the neutral ion fraction rises, photoabsorption heating becomes ef- 
ficient, and the gas enters a "Photoabsorption Plateau", consisting 
of WIM and WNM components, in which the temperature declines 
much more slowly. Finally the gas cools to the minimal computed 
temperature. See text. 

main nearly equal, leading to the slow thermal evolution 
in the plateau. However, departures from exact thermal 
equilibrium enable the gas to cool further, down to a 
temperature of 1650 K at tt-q x ^ ~ 5.4 x 10^^ cm~^ s. 

As the gas cools, and the heating radiation is ab- 
sorbed, the H and He preferentially absorb low energy 
photons, and the remaining ionizing radiation therefore 
becomes harder with increasing depth in the photoab- 
sorption plateau. Figure [3^ shows that after the UV 
photons are absorbed (at i ~ 5 x 10'^'* s), the heating 
rate decreases much more rapidly, and the cooling rate 
lags slightly behind the heating. This leads to the brief 
temperature dip at uq x t = 5.4 x 10^'^ s cm"'^ in Fig. 3a. 
After reaching a temperature minimum of 1650 K at this 
time, the gas heats up again to ~ 4500 K. After this, 
the cooling rate overcomes the heating rate again, and 
the gas then cools monotonically to 1000 K where we 
terminate the computation. 

The shock structure and time-scales that are shown in 
Figure [1] are in qualitative agreement with previous com- 
putations (Dopita & Sutherland 1996, Allen et al. 2008). 
For example, our results for the temperature profiles in a 
solar metallicity shock, are similar to those presented in 
Allen et al. (2008) for their "Dopita 2005" abundance set 
(middle panel of their Figure 7), where gas stays in the 
hot radiative phase for a few x 10^^ s, and then rapidly 
cools through the non-equilibrium cooling zone to enter 
the photoabsorption plateau. Differences between our re- 
sults and those of Allen et al. are likely due to differences 
in the assumed abundances, and in the strength of the 
magnetic field, B. 

3.2. Dependence on Gas Metallicity 

The time-scale over which the gas cools from the ini- 
tial shock temperature down to 1000 K, depends on the 
shock velocity that sets the initial energy content of the 
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shocked gas, and on the coohng efhciency of the post- 
shock gas. For = 5 x 10^ K, the coohng is initiaUy 
dominated by metal resonance hne coohng for Z > 0.1, 
and by bremsstrahlung emission for Z < 10^^. For 
Tg > lO'' K, bremsstrahlung cooling becomes important 
even for high-metallicity gas. Metal lines continue to 
dominate the cooling at lower temperature for Z > 0.1, 
whereas for Z < 0.1, hydrogen and helium dominate the 
cooling at T < 3 X 10^ K (see also GS07). 

The cooling efficiency therefore strongly depends on 
the metal content within the gas. For Z > 0.1, the cool- 
ing efficiency is roughly proportional to the gas metal- 
licity. However, at lower gas metallicity, as the metal 
contribution to the cooling rate above ~ 10^ K becomes 
negligible, the cooling efficiency approaches a limit set 
by the primordial helium abundance (see Boehringer & 
Hensler 1989; GS07). This can be seen in Figure [TJ For 
high gas metallicities the cooling time is proportional to 
Z, whereas the curves for Z = 10~^ and for Z = 10^"^ 
nearly overlap for temperatures above 10^ K. 

For Z > 1, even at ~ 10** K, permitted metal transi- 
tions dominate the cooling, and are much more efficient 
than Lya. For Z = 1 Lya provides ~ 15% of the cooling 
at its peak efficiency. For lower gas metallicities, Lya 
provides most of the cooling at T ~ 10^ K. For example, 
for Z = 0.1, Lya provides as much as 80% of the cool- 
ing. Therefore, the cooling times at ~ 10'' K are nearly 
independent of gas metallicity, except ior Z > 1. 

As the gas cools below ^ 10** K, metal fine structure 
emissions start to dominate the cooling, even for the low- 
est gas metallicities. In addition, the gas opacity and 
associated depth over which the heating radiation is ab- 
sorbed are functions of the gas metallicity. The cool- 
ing times below ^ 10* K are roughly proportional to 
gas metallicity for all Z. The overall total cooling times 
in isochoric flows is therefore sensitive to Z, as can be 
clearly seen in panel (a) of Figure [T] 

The different cooling times also set the degree of non- 
equilibrium ionization in the gas (GS07), as affected 
by the ratio of the cooling-time and the (metallicity- 
independent) recombination time. When the cooling is 
faster than recombination, departures from equilibrium 
may occur. Since cooling is faster for high gas metallic- 
ities, non-equilibrium effects, and the over-ionization in 
the gas, are more important for higher Z. These effects 
occur at temperatures between ^ lO^K and 10** K where 
the cooling is rapid. Below 10'* K, the heating becomes 
efficient, and the cooling times becomes long, so that the 
gas stays close to photoionization equilibrium. 

A second factor that is strongly affected by the gas 
metallicity is the radiation field emitted by the cool- 
ing gas. As discussed before, the emitted radiation is 
composed mainly of bremsstrahlung continuum, and of 
emission-lines and recombination continua. The total 
flux of radiation emitted by the cooling gas equals the 
input energy flux into the flow, 

F (X nokBTsVs cx uqv^ . (11) 

However, the relative contribution of metal emission lines 
to the total flux depends on Z. 

For high-metallicity g large fraction of the input 
energy is radiated as line emission. For low metallicity, 
the relative contribution of lines is small, and most of the 
initial energy flux is radiated as thermal bremsstrahlung. 



This can be clearly seen in Figure [2l For Z — 2 (the 
gray solid line), the lines-to-continuum contrast is large, 
producing the "UV-bump" discussed above. For Z = 
0.1 (dashed dark curve), lines are still important, but 
are much less pronounced than for Z — 2. For Z — 
10~^ (black solid curve) the spectrum is very smooth, 
and consists almost entirely of thermal bremsstrahlung. 

The differences in the spectral energy distribution of 
the emitted radiation field affect the thermal evolution 
of the cooling gas in the photoabsorption region, and the 
resulting ion distributions. We discuss the impact of the 
changing spectral energy distribution on the ion fractions 
in detail in § H) 

Additional features in the shock profiles are affected by 
the gas metallicity and associated cooling efficiency. The 
start of the photoabsorption plateau is associated with 
the increase in photoabsorption efficiency that occurs as 
the neutral hydrogen fraction becomes sufficiently large 
(of order 10^"^ for 5 x 10® K shocks). Higher-metallicity 
gas is more over-ionized, and therefore reaches this criti- 
cal neutral fraction at a lower temperature. The start of 
the WIM plateau therefore takes place at lower temper- 
ature for higher metallicity gas. For Z = 2, the plateau 
starts at 10* K, for Z = 1 at 1.7 x 10* K, and for Z < 0.1 
at 2.5 X 10* K. The depth of the WIM plateau is given 
by equation (|10p . Since the energy flux emitted by mod- 
els with the same initial shock temperature and magnetic 
field is similar, the depth of the WIM plateau is indepen- 
dent of gas metallicity. For our = 5 x 10® K shocks it 
is ~ 2.5 X 10^* s (for uq — 0.1 cm^"^ and strong-_B) for Z 
between 10^"^ and 2. 

For Z < 1, a, second warm ionized medium plateau 
is apparent in Figure [TJ This plateau has a tempera- 
ture T ^ 10* K. The depth over which the heating self- 
radiation is absorbed and the WNM plateau persists, de- 
pends on the gas opacity which is a function of Z. Lower 
metallicity gas produces longer WNM plateaus. 

In the final cooling of the gas, fine-structure transitions 
dominate the cooling. The intensity of fine-structure 
cooling is proportional to gas metallicity as discussed 
above. At high gas metallicities {Z > 1) fine-structure 
cooling is efficient, and starts to dominate the cooling at 
~ 10* K. At lower metallicities {Z < 0.1), fine-structure 
cooling starts to dominate at ~ 7500 K. 

3.3. Shock Temperature 

Panel (c) of Figure [1] shows the temperature profiles in 
a shock with an initial temperature = 5 x 10^ K, in 
the strong- _B limit. The overall characteristics of these 
temperature profiles are similar to those of the lower- 
velocity (Ts = 5 X 10® K) shocks discussed in § EH The 
gas is initially heated to Tg, then goes through a pro- 
longed hot radiative phase during which it emits most 
of its initial thermal energy as radiation. As T drops, 
cooling becomes more efficient, and a phase of rapid non- 
equilibrium cooling begins, bringing the gas to a temper- 
ature of a few xlO* K. At this point the neutral hydro- 
gen fraction becomes large enough that photoionization 
heating becomes significant, and a temperature plateau 
is formed. The WIM plateau ends when the neutral 
fraction is of order unity. A WNM plateau then fol- 
lows. The gas finally cools to our cut-off temperature, 
Tiow = 1000 K. 

Since the initial temperature is 10 times higher in a 
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5 X 10'' K shock than in a 5 x 10^ K shock, and the 
velocity is 10^/^ ~ 3 times higher, the total energy flux 
input is ^ 30 times higher. The overall cooling time is 
therefore longer, as the gas has to radiate more energy. 



Z = 2, "Strong-B" limit 




E (Ryd) 

Fig. 5. — Shock self-radiation as a function of shock temperature 
for the strong-S Umit, and for Z = 2. The dark curve is the mean 
intensity (erg s"'^ cm~^ Hz"'^) for Ts = 5x 10® K. The gray curve is 
the mean intensity for Ts = 5 X lO'^ K. The spectra are normaUzed 
so that their UV continuum-intensities overlap, to emphasize the 
different spectral energy distributions. 

The emitted radiation field is affected by the increased 
initial temperature. The higher temperature produces a 
more intense bremsstrahlung continuum with more ener- 
getic photons (the exponential bremsstrahlung cut-off oc- 
curring at higher photon frequencies). In addition, since 
at temperatures above ~ 2 x 10^ K bremsstrahlung is the 
dominant cooling process even at Z = 2, a large fraction 
of the initial energy is radiated as bremsstrahlung contin- 
uum rather than in lines. Lines still dominate the cooling 
at lower temperatures, but their total contribution to the 
integrated spectrum is smaller. This is shown in Figure[5] 
that compares the (normalized) spectral energy distribu- 
tions for 5 X 10^ K and 5 x 10^ K shocks. The hotter shock 
provides a flatter spectrum extending to higher energies, 
and its line-to-continuum contrast is smaller due to the 
larger fraction of radiation emitted at high temperatures 
dominated by bremsstrahlung emission. 

Because much of the initial energy in the 5 x 10^ K 
shock gas is radiated as bremsstrahlung, the cooling 
times above ~ lO"' K are much less sensitive to the gas 
metallicity than for 5 x 10^ K-shocks. In fact, the cool- 
ing time from 5 x 10^ K to ~ 10"* K increases only by a 
factor of 1.6 when the metallicity is reduced from Z = 1 
to Z = 0.1. 



the evolution below 10** K, while qualitatively similar to 
that of the previously discussed isochoric shocks, is com- 
pressed into a very short interval. This can be seen in 
Figure [TJs. Once the gas starts to cool, the decline down 
to a temperature of 1000 K is very rapid. 
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Fig. 6. — Shock structure versus time for Ts 



Z = 
tion, 



X 10 

: 5 X 10'' K, 

2 solar, for _B = 0. The plots focus on the final evolu- 
between 2.75 X IQi-* s and 2.8 X lO^"' s for no = 0.1 cm-3. 
(a) Temperature (K). (b) Electron fraction, (c) Neutral hydro- 
gen fraction, (d) Integrated H I column density (cm~'^) from the 
shock front, (e) Local cooling (dark) and heating (gray) rates per 
volume (erg s~^ cm""^). 



3.4. Magnetic Field 

Panel (b) of Figure [1] shows the temperature profiles for 
Ts = 5 X 10^ K, assuming B = 0. As discussed in § 
the hydrodynamics follow equations ([3]), and the evolu- 
tion is nearly isobaric, with = 4/3 Pq. The most 
important difference between shocks in which B = and 
the previously discussed strong-B limit, is that for B — 
the density increases as the gas cools. Since the cooling 
time is inversely proportional to the gas density, this im- 
plies a rapidly decreasing cooling time within the flow. 
The overall cooling times are therefore much shorter, and 



Figure [6] shows a "zoomed-in" snapshot of the final 
stages of the evolution for Z = 2. All times are shown 
assuming a post shock hydrogen density of 0.1 cm~'^. 
The insert in panel (a) shows the full temperature pro- 
file in the fiow. The panels focus on the final evolution, 
between t = 2.75 x lO^^^ s and 2.8 x lO^'' s. The final 
evolution shows similar features to those discussed for 
the strong- i? isochoric shocks. The gas cools rapidly to a 
point where photoabsorption becomes significant enough 
that the heating balances the cooling. It then enters a 
photoabsorption region in which the shock self-radiation 
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is gradually absorbed, and the temperature decline be- 
comes slower. This happens on much shorter time scales 
due to the increasing density (c.f. Figure 5 in Allen et 
al. 2008). 

Equation (fTT|) states that the total flux created by 
the cooling gas is proportional to the input energy flux, 
F cx nov^. Almost all of this flux (> 99.8%) is emit- 
ted before the plateau starts, and most of it (> 80%) is 
emitted within the hot radiative zone. The level of pho- 
toionization in the gas is determined by the ionization 
parameter, which is proportional to F/n. For strong- _B 
shocks, the density in the flow is constant, and the ioniza- 
tion parameter is therefore cx . For B = 0, the density 
increases as the gas cools, and the ionization parameter 
is therefore cx v^no/n, which is smaller by a factor no/n 
(or T/Ts). Photoionization is therefore much less im- 
portant in i? = models, and its contribution to the 
creation of intermediate- and high-ions is diminished. 



T =5x10^K, Z = 2 




10° 10^ 10' 

E (Ryd) 

Fig. 7. — Shock self radiation for Ts = 5 X 10^ K shocks and 
for Z = 2 solar metallicity gas. The solid curve shows the mean 
intensity (erg s~^ cm~^ Hz~^) in the strong-_B limit. The gray 
dashed curve shows the mean intensity for B = 0, divided by 5/3 
to remove the PdV contribution. As expected, the spectra are 
similar. 

A more subtle effect is related to the PdV work that 
appears in B = shock models, due to the gas com- 
pression. The work done on the cooling gas implies that 
the overall emitted radiation field is 5/3 times larger for 
B — Q shocks than for isochoric shocks (Edgar & Cheva- 
her 1986; GS07). The final total flux in strong-B flows 
is 3/2 riQkBTsVs, while for _B = it is 5/2 nak^TgVs- 
This is shown in Figure [71 However, the effect of the 
5/3-factor on the downstream ionization parameter is 
much smaller than the impact of the increasing density 
va B = Q shocks. 

Figure [8] compares the temperature profiles for strong- 
B and i? = models. Initially, due to the PdV work, the 
cooling in the isochoric, strong- i?, shock is faster than in 
the B = nearly isobaric flow. However, for B = Q 
the density and cooling rates quickly grow, and the final 
B = Q cooling is much more rapid. 

Figure [5] shows the hydrodynamic evolution of a i? = 0, 
T, = 5 X 10*^ K, Z = 2 fiow. Panel (a) shows the gas 
temperature. Panel (d) shows the gas pressure. As dis- 
cussed above (see equation [5]), the flow is nearly iso- 
baric, with a total pressure change of 4/3. This implies 




X t (cm ^ s) 



Fig. 8.— Temperature profiles for = 5 X 10'= K and Z = 0.1 
solar metallicity. The solid curve shows results in the strong-B 
limit. The dashed curve shows results for B = 0. 



that the gas density is roughly inversely proportional to 
the temperature, as can be seen in panel (b). The mass 
flux conservation equation implies that the velocity is in- 
versely proportional to the gas density as can be seen in 
panel (c). 

For B — Q, the pressure is dominated by the thermal 
gas pressure through the flow. We can therefore ver- 
ify mass, momentum and energy conservation, and we 
indeed find that they are conserved to better than 1%. 
Panel (e) shows the energy fiux components in the flow. 
The total energy flux, shown by the upper solid line, re- 
mains constant as it should. The thermal energy compo- 
nent {5Pv/2) is shown by the dark dashed line. Initially 
it dominates the energy flux, but as the gas cools, the 
thermal component gradually decreases. The bulk ki- 
netic energy {pv'^ /2) is shown by the dash-dotted line. 
Initially, this is only a small fraction of the total energy 
flux, and it decreases further as the gas cools and decel- 
erates. The gray line shows the integrated cooling ra- 
diation. As expected, this component grows with depth 
into the flow, and finally reaches a value equal to the to- 
tal initial energy flux when all internal thermal energy is 
lost. 

The impact of gas metallicity on B = isobaric shocks 
is similar to that which occurs in strong-i? isochoric 
shocks discussed in § 13.21 Again, as the gas metallic- 
ity increases, the metal-line cooling efficiency increases, 
and cooling becomes more rapid. This is clearly seen in 
Figure [T]d, where the initial cooling time is proportional 
to gas metallicity, except at Z < 10~^ where the metals 
contribution to the cooling is negligible. 

The cooling below ~ 10^ K is very rapid due to the 
increased gas density. While the cooling times below 
^ 10'' K are proportional to gas metallicity, this cannot 
be seen in Figure [TJd since the total cooling times below 
'--^ 10^ K are so much shorter than the cooling times from 
5 X 10® K to 10^ K. The low-temperature cooling time 
only becomes long enough to be seen in Figure [TJd for 
Z < 10~^, for which the cooling time during the plateau 
phase is of order the cooling time to reach the plateau. 

3.5. Gas Density 
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Fig. 9.— Dynamics in a 5 X 10*^ K, Z = 2, B = shock, (a) Tem- 
perature (K). (b) Density profile, (c) Velocity profile, (d) Pressure 
profile, (e) Energy flux (erg s~^ cm~^). The dark solid curve 
shows the total energy fiux in the flow, which is conserved as it 
should. The dashed line shows the thermal energy fiux {^Pv), 
which is initially the dominant component of the energy flux, but 
drops as the gas cools. The dash-dotted line shows the bulk kinetic 
energy fiux (^pi)^). The gray solid curve is the integrated cooling 
radiation, which is accumulated in the flow as the gas cools, until 
it reaches the total input energy fiux when all thermal energy is 
lost. 



The results discussed above for T versus tiq x i are inde- 
pendent of density for strong-B isochoric models. When 
B = 0, the gas is compressed as it cools, and the thermal 
evolution below ^ 10^ K may depend on the initial den- 
sity. In the high-density compressed gas, the cooling ef- 
ficiencies may be suppressed by collisional de-excitations 
of the cooling transitions if the densities become suffi- 
ciently high. The temperature at which such collisional 
quenching occurs depends on the shock velocity as well 
as on the initial post-shock density, hq. 

Figure [To] shows the temperature as a function of no x t 
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Fig. 10. — Temperature profiles for different post-shock hydro- 
gen densities, for a 5 X 10*' K shock and B = 0. The plot focuses on 
the final evolution. Higher density models cool less efficiently be- 
cause collisional quenching occurs sooner as the gas is compressed. 

in Ts = 5 X 10® K shocks, assuming post-shock hydrogen 
densities of 10~^ cm~^ (dashed line), 0.1 cm~^ (solid 
line), and 100 cni^^ (dash-dotted line). For these ini- 
tial densities, the gas density never becomes high enough 
for significant collisional de-excitation of the permitted 
transitions. However, once the temperature drops below 
~ lO'' K, some of the fine-structure cooling is quenched, 
and the thermal evolution is altered. As shown by Fig- 
ure [ini for higher initial post-shock densities this colli- 
sional quenching occurs sooner, and the effective cooling 
times, no x t, are therefore longer for higher no (c.f. Fig- 
ure 6 in Allen et al. 2008, in which the cooling times from 
Tg to the photoabsorption plateau scale as n~^, whereas 
the cooling below ~ 10"* K depends more strongly on n). 

3.6. Steady-State Conditions 

The results presented in this paper rely on the assump- 
tion that the shocks have reached a steady-state struc- 
ture. Attaining steady-state requires that the shock ex- 
ists for a time-scale that is longer than the cooling time. 
Otherwise, the shocked gas only cools partially, and the 
shock structure, self-radiation, and integrated column 
densities are time-dependent. 

In Table [2] we list the total cooling times (yr), from 
down to 1000 K, assuming a post-shock hydrogen den- 
sity, no = 1 cm"'^. We also list the associated length- 
scales (kpc), and total hydrogen column densities, A'n.tot 
(cm^^). These are then the hydrogen columns and length 
scales required for the formation of steady-state shocks. 
For example. Table [2] imphes that for a Tg = 5 x 10® K 
shock in solar metallicity gas, the time required to reach 
steady-state is 6.4 x 10® yrs, in the strong-i? isochoric 
limit. The cooling length is then 1 kpc, and the asso- 
ciated hydrogen column density is 3 x 10^^ cm^'^. The 
time- and size-scales are generally proportional to Uq^, 
while the total hydrogen columns are generally indepen- 
dent of no (but see § 13. 5p . 

4. ION FRACTIONS 

We have computed the ionization states of H, He, C, 
N, O, Ne, Mg, Si, S, and Fe in the post-shock cooling 
layers. When the photoionized precursor-gas enters the 
shock, its ionization state rapidly reaches CIE at a tem- 
perature very close to T,. As the hot gas flows away from 
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TABLE 2 

Shock Parameters (for no = 1 cm-^) 







Ts 


= 5 X 10 


6 K 


Ts 


= 5 X 10 


' K 




Metallicity 


Time 


Size 


A'H,tot 


Time 


Size 


A^H.tot 






(yr) 


(kpc) 


{cm-2) 


(yr) 


(kpc) 


(cm"2) 


strong-_B (isochoric) 














Z = 


2 


2.9 X 10° 


0.45 


1.4 X 102i 


9.2 X 10' 


45 


1.4 X 1023 


Z = 


1 


6.4 X 10^ 


1.0 


3.0 X 10^1 


1.4 X 10® 


71 


2.2 X 1023 


Z = 


10-1 


5.7 X 10'' 


8.8 


2.7 X 1022 


5.2 X 10* 


250 


7.8 X 1023 


Z = 


10-2 


3.3 X 10** 


50 


1.6 X 1023 


1.5 X 10^ 


750 


2.3 X 1024 


Z = 


10-3 


2.1 X 10^ 


330 


1.0 X 1024 


4.1 X 10^ 


2000 


6.3 X 1024 


B = 


("isobaric") 














Z = 


2 


8.9 X 10^ 


0.10 


4.3 X 1020 


2.2 X 10^ 


7.0 


3.3 X 1022 


Z = 


1 


1.6 X lO'' 


0.18 


7.8 X 1020 


2.8 X 10^ 


8.6 


4.2 X 1022 


Z = 


10-1 


6.1 X 10^ 


0.63 


2.9 X 1021 


3.8 X lO'' 


11.0 


5.7 X 1022 


Z = 


10-2 


1.0 X 10'' 


0.87 


4.8 X 102i 


4.1 X 10^ 


11.4 


6.3 X 1022 


Z = 


10-3 


1.6 X 10'' 


0.91 


7.7 X 102i 


5.3 X 10^ 


11.5 


8.0 X 1022 



Note. — Cooling-times (yr), size-scales (kpc), and total hydrogen column densities 
(cm— 2) for post-shock gas cooling from Tg to 1000 K, assuming a post-shock hydrogen 
density, ng = 1 cm" 3. 



the shock front, it recombines, cools, and radiates away 
its thermal energy. This radiation is later absorbed by 
the cooler gas further downstream, providing a source of 
heating and photoionization. Wc follow the time depen- 
dent ion fractions in the flow, taking into account pho- 
toionization by the shock self-radiation. As we discuss 
below, the photoionizing radiation significantly affects 
the ion fractions in the gas. 

When the cooling time becomes short compared to the 
recombination time, departures from equilibrium may 
occur, keeping the gas over-ionized compared to ion- 
ization equilibrium at the local conditions (as specified 
by the mean intensity of the photoionizing radiation, 
gas density, and temperature). We consider the non- 
equilibrium ionization states as a function of the time de- 
pendent temperature and mean radiation intensity. We 
present results for gas cooling behind shocks with initial 
post-shock temperatures 5 x 10^ K and 5 x 10^ K, with 
metaUicities Z = 10~^, 10~^, 10~^, 1, and 2 times the 
solar metal abundances, in the B = (nearly isobaric) 
and "strong-B" (isochoric) limits. 

Table |3] lists the ion fractions as a function of time 
and temperature for the various models that we consider, 
as outlined in Table S] All the times given in Table [3] 
were computed assuming a post-shock hydrogen density 
nufl = 0.1 cm"^. As we discussed in §[3l the shock struc- 
tures, as functions of no x t, are independent of the gas 
density in the strong-B limit, and are density-dependent 
only at low temperatures (< 10^ K) when B — 0. 

As a detailed example, in Figure[TT]wc show the carbon 
ion fractions for = 5 x 10^ K, in the strong-_B limit, 
for the various values of Z. The second panel (Fig [TTb) 
shows the results for Z = 1. We first focus on these 
results, and then consider how the results depend on Z. 

As the photoionizcd precursor gas enters the shock, the 
carbon ion fractions approach CIE at the shock temper- 
ature. This initial rise is very rapid, and cannot be seen 
in Figure [TlJ Figure [12] focuses on this initial evolution. 
The abundant carbon ions in the precursor gas are C^"*" 
and C^"''. As can be seen from Figure [121 these ionization 
states are quickly replaced by C^+, which is the abun- 



dant ion at = 5 X 10^ K (see GS07). By the time 
the temperature drops to 4.97 x 10^ K, the carbon ion 
fractions reach CIE. Later on, as carbon recombines in 
the cooling flow, the lower ionization states reappear at 
lower temperatures, as can be seen in the second panel 
of Figure [TT] For example, is again the most abun- 
dant carbon ion between T - 10^ K, and T - 2 x 10"* K. 
Many ionic species show such a double peak abundance 
pattern, where the abundance peaks first immediately 
after passing through the shock front, and then a second 
time as the gas recombines further downstream. 

The time scales of the two peaks are very different. 
For the example shown is Figure [T^ the width of the 
first peak in the C"'"*" abundance is of order 10^^ s, while 
the second peak lasts of order 10^'' s. The total contribu- 
tion of the first abundance peak to the total ionic column 
density is therefore very small. However, since the gas is 
underionized as it goes through this first peak, the hot 
ambient electrons efficiently excite the transitions of the 
underionized ions (McCray 1987), creating an enhanced 
emissivity. This enhanced emissivity also implies that 
the cooling rates during this adjustment phase are en- 
hanced. In lower velocity shocks, the first peaks can con- 
tribute significantly to the integrated column denisties 
(Krolik & Raymond 1985). 

Figure [Tl] shows that in the hot radiative zone, the 
gas remains close to CIE, with C^"*" being the most dom- 
inant carbon ion. During the non-equilibrium cooling 
phase (T < 10^ K), the gas cools and recombines, and 
the most abundant ionization state drops to C^^ and 
later to C^^. The ionization state stays higher than at 
CIE due to photoionization by the shock self-radiation 
as we discuss below. The photoabsorption plateau starts 
at T '--^ 2 X 10^ K, and as the shock self radiation is ab- 
sorbed within the plateau, the dominant carbon ion grad- 
ually drops from C'^"'" to C^"*" , and finally to C"*" . For the 
parameters considered here, C"*" remains the dominant 
species until the gas reaches our termination temperature 
21ow = 1000 K. Our results are in qualitative agreement 
with those of Allen et al. (2008; see their Figure 9), in 
which C^+ is the most dominant carbon ion in the hot ra- 
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TABLE 3 

Ion Fractions in a 5 x 10^ K, 2 times solar, strong- _B-shock 



Time 


Temperature 


HO/H 


H+/H 


He'' /He 


s 


K 








0.00 


5.000 X 10^ 


1.63 X 10-* 


0.00 


4.93x10-3 ... 


7.25 X 10^ 


4.998 X 10^ 


1.52 X 10-^5 


0.00 


3.59x10-3 ... 


1.57 X 10^ 


4.997 X 10^ 


1.01 X 10-*^ 


0.00 


2.48x10-3 ... 



Note. — The complete version of this table is in the electronic edition 
of the Journal. The printed edition contains only a sample. The full 
table lists ion fractions for the B = (isobaric) and strong-B (isochoric) 
magnetic field limits, for shock temperatures of 5 X 10® K and 5 X 10^ K, 
and for Z = 10-3, IQ-^, IQ-^, 1, and 2 times solar metallicity gas (for 
a guide, see Table |4jl. The times in the first column are fo an assumed 
post-shock hydrogen density of 0.1 cm— 3. 



TABLE 4 

Ionization and Cooling Tables 



Data 








Ion Fractions 


Cooling 


Post-Shock 
Columns 


Precursor 
Columns 


5 X 10'^ 


K 


strong- _B 


Z = 2 


3A 


5A 


6A 


8A 








Z = 1 


3B 


5B 


6B 


8B 








Z = 10-1 


3C 


5C 


6C 


8C 








Z = 10-2 


3D 


5D 


6D 


8D 








Z = 10-3 


3E 


5E 


6E 


8E 






B = 


Z = 2 


3F 


5F 


6F 


8F 








Z = 1 


3G 


5G 


6G 


8G 








Z = 10-1 


3H 


5H 


6H 


8H 








Z = 10-2 


31 


51 


61 


81 








Z = 10-3 


3J 


5J 


6J 


8J 


5 X 10^ 


K 


strong-B 


Z = 2 


3K 


5K 


6K 


8K 








Z = 1 


3L 


5L 


6L 


8L 








Z = 10-1 


3M 


5M 


6M 


8M 








Z = 10-2 


3N 


5N 


6N 


8N 








Z = 10-3 


30 


50 


60 


80 






B = 


Z = 2 


3P 


5P 


6P 


8P 








Z = 1 


3Q 


5Q 


6Q 


8Q 








Z = 10-1 


3R 


5R 


6R 


8R 








Z = 10-2 


3S 


5S 


6S 


8S 








Z = 10-3 


3T 


5T 


6T 


8T 



diative phase, C^^-C"^^ dominate in the non-equilibrium 
cooling zone, and in the plateau C^"*" recombines to from 
C+, which remains dominant down to 1000 K. 

4.1. Photoionization 

The radiation emitted by the hot post-shock gas has 
a profound effect on the ion fractions in the cooling gas. 
To illustrate the importance of photoionization by the 
gas self-radiation on the evolution and ion fractions, we 
have computed a comparison-model in which the radia- 
tion field at any point in the flow is artificially set to zero, 
so that there is no photoionization and associated heat- 
ing anywhere in the flow. The results (for Tg — 5x 10^ K, 
strong-B, and Z = 1) with and without photoionization 
are presented in Figure 1131 

In the hottest parts of the flow where T > 10^ K, 
the ionizing radiation has a negligible effect on the ion 
fractions, because the typical ionization potentials are 
too high to be efficiently affected by the radiation. At 
lower temperatures, photoionization becomes important. 
For example, in the post-shock gas the ion frac- 
tion remains higher than 10% down to a temperature of 
~ 2 X 10* K, whereas without photoionization it drops 



below 10% at T ~ 10^ K. As expected, photoionization 
becomes dominant at lower temperatures, maintaining 
as an abundant species below Iff* K, and 0+ the 
most abundant species even at T ~ 10^ K. Without ra- 
diation, the ionization level is much lower, and neutral 
carbon is the most dominant species below 5000 K. Pho- 
toionization thus strongly affects the resulting integrated 
column densities in the cooling layer. The column den- 
sity of the high-ion is enhanced by a factor of ~ 1.5 
due to photoionization. The column densities of the mid- 
and low-ions C^+, C^+, and C+ are enhanced by factors 
of more than 100. 

4.2. Shock Temperature 

The shock self-radiation depends on the shock veloc- 
ity and associated initial post-shock temperature Tg. As 
discussed in § 13.21 (equation [H]), the mean intensity of 
the radiation field is proportional to T^ '^. In addition 
to the overall intensity dependence, the spectral-energy 
distribution hardens with increasing Ts . Figure [5] shows 
the spectral energy distributions for Tg equal to 5 x 10^ K 
and 5 x 10^ K. Even high ions, which are only collision- 
ally ionized in a 5 x 10^ K shock, are photoionized by 
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Fig. 11. — Carbon ion fractions versus gas temperature in a 
Ta = 5 X 10® K shock in the strong-S limit. The various panels 
show results for gas metallicitics ranging from 10^^ to 2 times 
solar. 
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Fig. 12. — Carbon ion fraction versus gas temperature just after 
crossing the shock front, for Ts = 5 X 10" K, strong-S, and Z = 1 
gas. The plot shows how the initial ionization states quickly adjust 
to CIE at the initial shock temperature. 

the harder photon emitted in a 5 x 10^ K shock. The 
gas in the hotter shock is more highly ionized, both due 
to the higher intensity, and the harder spectral shape of 
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Fig. 13. — Carbon ion fractions versus gas temperature in a 
Ts = 5 X 10® K strong-B shock for Z = 1. The upper panel shows a 
non-physical model in which photoionization is not included in the 
computation. The lower panel shows results when photoionization 
is included. 

the shock radiation. 

In Figure [TJ] we show the carbon ion fractions as a 
function of gas temperature. The upper panel is for — 
5 X 10'' K, and the lower panel is for Tg = 5 x 10^ K. It 
is clear that the gas is more highly ionized for the hotter 
shock. This is evident even at very high temperatures, 
> 10^ K, where the more energetic photons created by 
the hotter bremsstrahlung continuum of the 5 x 10^ K- 
shock efficiently ionize to C^"''. 
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Fig. 14. — Carbon ion fraction versus gas temperature in a solar 
metallicity gas, for a strong- iJ (isochoric) shock. The upper panel 
shows results for a shock temperature of 5 X 10^ K, and the lower 
panel shows results for 5 X 10® K. 

Some of the differences between the ion distributions 
shown in Figure [14] for the two shock velocities, result 
from the fact that the various "zones" occur at different 
gas temperatures depending on Tg. For example, the 
ranges of temperatures at which the gas is in the WIM 
plateau (see § [3]) are marked in Figure [U] by the gray 
shaded zones. The WIM plateau starts at a temperature 
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of ~ 5.5 X 10"* K for = 5 X 10'' K, and at a temperature 
of - 1.6 X 10-* K for = 5 X 10^ K. As the radiation 
is absorbed in the plateau, the carbon ionization states 
drops to C'^+, then C^+, and eventually to C"*" towards 
the end of the WIM plateau where hydrogen and helium 
become neutral. While this recombination process takes 
place within the WIM plateau for both values of Tg , the 
ion distributions versus temperature, XiiT), are different. 

In the absence of photoionization, the ion fractions, for 
ions that are produced collisionally at temperatures less 
than Ts, are independent of the shock velocity. Thus, 
without photoionization the ion fractions versus temper- 
ature shown in Fig. 17 would be identical for the two 
shock temperatures. The integrated column densities 
through the flow are then proportional to the shock ve- 
locity (e.g. Heckman et al. 2002). However, as is clearly 
seen from Figures [T51 and [HI photoionization plays a ma- 
jor role in setting the ion fractions in the flow, and the 
ion distributions are affected by the metallicity, shock 
temperature, and magnetic field, through their control 
of the ionization parameter. We present detailed results 
for the ionic column densities in § [6l 

4.3. Magnetic Field 

One of the parameters that determines the level of pho- 
toionization in the post-shock gas is the strength of the 
magnetic field. As discussed in § 13.41 the magnetic field 
strongly affects the ionization parameter in the down- 
stream gas, due to the compression that takes place when 
B is small. We therefore expect that at low tempera- 
tures, after significant compression has taken place, mod- 
els with B ~ Q will be much less ionized than strong-S 
isochoric models. 

Figure [12] shows the carbon ion fractions in the iso- 
baric 5 = (upper panel) , and isochoric strong-B (lower 
panel) limits. The lower ionization parameter in the 
downstream gas when B = Q results in lower ionization 
states at all temperatures. For example, at a tempera- 
ture of 30, 000 K, > 98% of the carbon is in the state 
of C'^+-C^+ in the limit of strong-i3, whereas for i? = 
more than half of the carbon is in the state of C°-C^+. 
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Fig. 15. — Carbon ion fractions versus gas temperature in a 
Ts = 5 X lO** K shock for Z = 0. The upper panel is for S = 0, 
and the lower panel is for strong- _B. 



For _B = 0, the temperature versus time in the post- 
shock cooling layers depends on the density, at temper- 
atures lower than ~ 10^ K (see § 13. 5p where collisional 
de-excitations of fine-structure transitions become signif- 
icant. The physical conditions below T ^ lO'* K, there- 
fore depend on the density of the shocked medium. The 
computations presented here assume a post-shock hydro- 
gen density of 0.1 cm~'^. 

4.4. Departures from Equilibrium Ionization 

In a photoionized gas with a high ionization parameter, 
the ion fractions depend only weakly on the gas tempera- 
ture, as opposed to purely collisionally ionized gas. Since 
the gas does not have to significantly adjust its ionization 
state to the time-dependent temperature during cooling, 
departures from equilibrium are expected to be smaller 
than for pure radiative cooling. 

To demonstrate how departures from equilibrium affect 
the ion fractions, we compare our time-dependent models 
to a computation in which local equilibrium ionization is 
imposed everywhere, given the local mean intensities and 
temperatures that obtain for 5 x 10^ K, Z = 1 shocks. In 
the comparison calculations we use the mean intensities 
and temperatures derived from the full time-dependent 
computations. (These mean intensities arc different from 
those that would be derived in a "self-consistent" equilib- 
rium model.) Figure [TB] shows the time-dependent (dark) 
and equilibrium (gray) carbon ion fractions as a function 
of temperature. The upper panel is for a strong-_B iso- 
choric model, and the lower panel is for 5 = 0. 

For strong-B shocks, the upper of Figure [TCj shows 
that departures from photoionization equilibrium occur 
at temperatures between ~ 2 x 10^ K and ~ 10^ K, but 
the ion fractions differ by < 25%. These differences are 
significantly smaller than the differences in the absence of 
photoionization (GS07). In strong- £? shocks, photoion- 
ization increases the abundances of high-ions at low tem- 
peratures, beyond the enhancement of the "collisional" 
recombination lag (see Figure [T3|). The time-dependent 
ionization states then remain close to photoionization 
equilibrium (to within 25%) even in the non-equilibrium 
cooling zone. 

The lower panel shows that when B — 0, departures 
from equilibrium ionization are larger, due to the smaller 
effect of photoionization resulting from the gas compres- 
sion (see 14. 3p . For example, the non-equilibrium abun- 
dance of 0"*+ is greater than 0.05 down to ^ 2.2 x 10"* K, 
whereas the equilibrium abundance vanishes below ~ 
3.3 X 10"* K. Departures from equilibrium ionization tend 
to keep the gas at any temperature over-ionized, as re- 
combination lags behind cooling. The contribution of 
photoionization increases the equilibrium ion fractions 
relative to the case of pure radiative coohng (GS07), 
especially at T < lO^K, where the shock self-radiation 
is energetic enough to efficiently ionize the abundant 
species, but the compression is still not too large to main- 
tain a significant ionization parameter. At lower temper- 
atures, the non-equilibrium ion-fractions are due entirely 
to the recombination lags in the collisionally ionized gas. 
Below ^ 10* K, efficient heating increases the net cool- 
ing times significantly, and the gas approaches photoion- 
ization equilibrium as can be seen by the near overlap 
between the dark and gray curves. 
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Fig. 16. — Carbon ion fractions versus gas temperature in Ts = 5 X 10® K shocks, for a solar metallicity gas. The upper panel is for the 
strong- limit, and the lower panel is for B = 0. The gray curves show the local photoionization equilibrium ion fractions. 



4.5. Gas Metallicity 

Our results for the non-equilibrium ion fractions, 
s^T), for metallicities Z equal to lO'^, 10"^ 10-\ 1 
and 2, are presented in Table [31 Figure [TT] shows, as an 
example, the carbon ion fractions for the different val- 
ues of Z. The assumed metallicity affects the ion frac- 
tions in several ways. First, the cooling times depend on 
the metal abundance. Higher Z leads to enhanced metal 
lines cooling, and therefore shorter cooling times. Depar- 
tures from equilibrium and recombination lags are there- 
fore larger for higher metal abundances. Higher metal- 
licity gas will tend to be more over-ionized. 

Second, the spectral shape of the shock self radiation 
depends on the gas metallicity. Since metal line emis- 
sion is enhanced for higher Z, a larger fraction of the 
initial energy is radiated via line emission, and a smaller 
fraction as bremsstrahlung continuum. High-metallicity 
shocks therefore produce more photons with energies be- 
tween 1 and ^ 5 Rydbergs, and fewer photons with 
^ ^ 6 Rydbergs. At low Z the shock self-radiation is 
harder. The changing spectral energy distribution as 
a function of gas metallicity affects the photoionization 
rates, and therefore the ion fractions as a function of 
temperature. 

As examples illustrating the various effects, in Fig- 
ure [IT^-c we display the 0"^+, 0^+, and Ne''+ distribu- 
tions for the different values of Z for a 5 x 10^ K shock 
in the strong-B limit. For comparison, in panels (d)-(f) 
we shows the local equilibrium ion fractions, given the 
local mean intensities and temperatures that obtain in 
the shock. 

Panel (a) shows that carbon is more over-ionized for 
higher metallicity gas, and that the ion fraction 

peaks at lower temperatures for higher Z. For exam- 
ple, for Z = 2 the distribution peaks at T ~ 10^ K. 
For smaller Z , the ion fraction peaks at higher temper- 
atures. For Z = 10"3, it peaks at T ~ 2 x lO'* K. The 
sharp decline in the C'^"'' fraction at ~ 10'' K, is related to 
the onset of the photoabsorption plateau, which occurs 
at a metallicity-dependent temperature (see § I3.2[) . Pan- 



els (a) and (d) show that the C'^^ fraction remains close 
to its equilibrium distribution for any value of Z, and for 
the entire temperature range over which it is abundant 
(- 10^^-2 X 10^ K). 

O^^ shows similar behavior. The higher metallicity 
shocks are more over-ionized, and the 0^+ fraction per- 
sists to lower temperatures where the photoabsorption 
plateau occurs. A comparison of panels (b) and (e) shows 
that departures from equilibrium ionization occur only 
for Z > 1, and are limited to T ~ 3 x 10^ K where the 
role of photoionization is still minor. 

For Ne^+, the behavior is a bit more complicated. At 
high gas metallicities {Z > 1), the recombination lags 
enhance the Ne^"*" abundances above the enhancements 
due to photoionization, for temperatures between 4 x 10^ 
and 2 x 10^ K. Because the spectral energy distributions 
for Z = 2 and for Z — I are similar, the equilibrium dis- 
tributions are identical. They are also narrower than the 
non-equilibrium distributions. However, for lower values 
of Z (< 0.1), the varying spectral energy distributions 
affects the Ne''+ fractions. Figure [2] shows that there are 
more photons capable of ionizing Ne^+ (with an ioniza- 
tion threshold of 15.2 Ryd) for low Z, than for Z = 2. 
Indeed, for decreasing Z the Ne^"*" ion distribution be- 
comes broader again, due to enhanced photoionization 
by the shock radiation. This can be seen in panel (f). A 
comparison with panel (c) then shows that for Z < 0.1 
the Ne^"*'-fraction approaches the equilibrium distribu- 
tion controlled by photoionization. 

At low gas metallicities (< 10^^), the total contribu- 
tions of the metals to the gas cooling and to the emis- 
sivity becomes negligible (GS07). The cooling rates and 
spectral energy distributions therefore become indepen- 
dent of Z. This can be clearly seen in Figure I17[ where 
the ion distributions for Z — 10~^ and for Z = 10~^ are 
similar for C^+, 0^+, and Ne''+. 

4.6. Auger Effects 

We have included multi-electron Auger ionizations us- 
ing cross section and yields from Kaastra & Mewe (1993). 
Auger ionization is a potential source of multiply charged 
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Fig. 17. — C'^'*", O^"'', and Ne'^^ ion fraction versus temperature for gas metallicities between 10"'^ and 2 times solar, in a 5 X 10^ K 
shock in the limit of strong-B (panels [a]-[c]). Panels (d)-(f) show, for comparison, the local equilibrium ion fractions, given the local mean 
intensities and temperatures that obtain in the shock. 



metal species in largely neutral (and cool) hydrogen gas, 
due to the large penetration depths of high-energy pho- 
tons. However, we find that rapid charge-transfer neu- 
tralization with atomic hydrogen quickly removes such 
ions. The metal abundances are therefore not signifi- 
cantly affected by the Auger processes (cf. Dopita and 
Sutherland 1996). However, the charge transfer reactions 
do lead to a slight rise in the ionized-hydrogen (proton) 
densities. We demonstrate this effect in Figure [T8l for a 
strong-B, Z = 1, 5x 10^ K shock, in which intense X-rays 
are produced. The dark and gray solid lines show the ion- 
ized hydrogen fractions with and without the inclusion 
of Auger processes. Major contributors to the increase 
in the ionized hydrogen fraction are charge-transfer reac- 
tions with Fe and Si ions. The dark dashed curve shows 
the H+ abundances if Auger ionizations are included for 
iron ions only. 




T (K) 



Fig. 18. — H+ fraction as a function of temperature for a 
strong-B, 5 X 10^ K shock at Z = 1. The dark solid line shows the 
H+ fraction when multi-electron Auger ionizations are included. 
The gray solid line shows a model in which Auger ionizations are 
excluded. The dashed curve shows the H+ ion fraction when Auger 
ionizations are only included for iron ions. 



5. HEATING AND COOLING 

We have carried out computations of the heating and 
cooling efficiencies in the post-shock cooling layers for 
shock temperatures of 5 x 10^ K and 5 x 10'' K, for 
gas metallicities Z between lO^'^ and 2, for B = and 
strong- i? dynamics. Our results for the heating rates and 
cooling efficiencies are listed in Table [5l as outlined by 



TABLE 5 

Cooling and Heating in a 5 x 10^ K, 2 times solar, 

STRONG- B- SHOCK 



Time 


Temperature 


A 


T 






(T) 


(erg s"-"- cm'') 


(erg s-i 


) 


0.00 


5.000 X 10^ 


2.72 X 10-21 


2.98 X 10" 


-27 


7.25 X 10** 


4.998 X 10^ 


2.34 X 10-21 


4.38 X 10" 


-28 


1.57 X 10« 


4.997 X 10** 


2.10 X 10-21 


1.79 X 10" 


-28 



Note. — The complete version of this table is in the elec- 
tronic edition of the Journal. The printed edition contains 
only a sample. The full table lists the cooling and heat- 
ing for strong-B and for B = 0, for shock temperatures of 
5 X 10^ K and 5 X 10^ K, and for Z = 10'^, 10'^, 10'^, 
1, and 2 (for a guide, see Table (Jjl. The times in the first 
column are given assuming a post-shock hydrogen density 
of 0.1 cm— ^. A and T are defined such that the cooling rate 
per volume is neUnA, and the heating rate per volume is 
ntotT. 

Table |4l We use the self-consistent non-equilibrium ion 
fractions, Xi, obtained in §[4] to compute A{T, Xi, Z) and 
T{xi, Z, Jy) for all temperatures in the downstream gas. 
There is no external source of photoionization and heat- 
ing, and all of the heating radiation is produced by the 
shock itself. In all of our shock models, photo-heating is 
negligible above temperatures of a few xlO^ K and the 
gas undergoes pure radiative cooling. However, as dis- 
cussed in § m photoionization of the metals can become 
a very important process even in the hotter components, 
and this can alter the radiative cooling efficiencies. 

For Z > 0.1, the radiative cooling between a few 
xlO^ K and ~ 10'' K is dominated by electron impact 
excitation of resonance line transitions of metal ions. 
Above ^ 10^ K, bremsstrahlung cooling dominates for 
all gas metallicities. At lower temperatures (< 10'* K), 
Lyof cooling and fine-structure lines dominate the cool- 
ing, with relative contributions depending on T and Z, 
as we discuss below (Boehringer & Hensler 1989; Suther- 
land & Dopita 1996; GS07). 

As an example, in Figure [T9l we show the cooling and 
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heating rates per unit volume, UeUH^ and riT, as a func- 
tion of time (in panel [a]), and as a function of tempera- 
ture (in panel[ b]) for a strong-B, Tg = 5 x 10^ K shock 
at Z = 1. In panel (a), we have assumed a post-shock 
hydrogen density of 0.1 cm^'^ for the cooling time-scales. 
The dominant cooling and heating elements are indicated 
near the curves. 

The gas starts out in the "hot radiative" phase. In 
this model it remains in this phase for ~ 3.5 x lO^"* s, 
and cools gradually, via metal resonance transitions and 
bremsstrahlung emission, down to a temperature of ~ 
2.5 X 10^ K. It then enters the "non-equilibrium cool- 
ing" phase, where it rapidly cools to a temperature of 
~ 10^ K via very efficient metal line cooling. This cool- 
ing is due mainly to iron, neon, and oxygen ions, as in- 
dicated in Figure [T9b . Throughout the "hot radiative" 
and "non-equilibrium cooling" phases, heating is negli- 
gible, and the gas cools radiatively from the initial hot 
shock temperature down to ^ 10"* K. However, since pho- 
toionization has a major influence on the ion fractions, 
the non-equilibrium cooling rates differ significantly from 
those presented in GS07 in which no sources of ionizing 
radiation were considered. We discuss this point further 
below. 

As discussed in § [3l once the neutral hydrogen frac- 
tion becomes large enough, > 10~^, to allow for efficient 
photoabsorption, the gas enters the "WIM photoabsorp- 
tion plateau" . In this stage, heating is mainly due to 
H and IIe+ photoabsorption. As previously discussed, 
for Z = 1, Lya provides only 15% of the cooling at its 
peak efficiency. The rest of the cooling is provided by 
metal resonance lines for T > 2 x 10** K, and by fine- 
structure (FS) and resonance line transitions at lower 
temperatures. For lower gas metallicities, Lya consti- 
tutes a larger fraction of the cooling at T ~ 10^ K. Be- 
low < 8000 K, Lya cooling vanishes and the cooling is 
dominated entirely by fine-structure transitions. 

Heating is first dominated by hydrogen photoabsorp- 
tion. As the radiation is increasingly absorbed and the 
gas cools to ~ 6500 K, helium photoabsorption dom- 
inates the heating. Finally, below ~ 2000 K, and as 
the ionizing radiation becomes even harder, oxygen and 
other metals become significant absorbers. In the WIM 
and WNM plateaus the gas remains close to thermal 
equilibrium, as is evident by the fact that the heating 
and cooling curves nearly overlap. The delicate balance 
between heating and cooling that depends on the tem- 
perature, electron density and ion fractions, determines 
the temperatures profiles that were discussed in § [3l 

The way in which the cooling efficiency. A, depends 
on the gas metallicity determines the column-density ra- 
tios in the post-shock cooling layers, as we discuss in 
§ |6] below. For T > 10^ K, cooling is dominated by 
bremsstrahlung emission for all gas metallicities, and the 
cooling efficiency is therefore only weakly dependent on 
Z. 

For 2 X 10* K ^ r < 10^ K, resonance metal line tran- 
sitions are efficiently excited. For Z > 0.1, the dominant 
cooling process is metal line cooling, and A cx XjZ. For 
Z < 0.01, the contribution of metals to the gas cooling is 
negligible. The cooling is dominated by bremsstrahlung 
emission and by hydrogen and helium Lya, and A is 
therefore independent of Z. 

At T ^ 10* K, a mixture of metal fine structure and 



TABLE 6 

The Dependence of A{xi,T) on Z 



Temperature 


Cooling efficiency, A{xi 




K 


Z < 10-2 


10-2 < Z <1 


Z > 1 


< 7000 


ocZ (1) 


oc Z (1) 


oc Z (1) 


~ lO* 


ocl (2) 


cxl (2) 


ocZW 


2 X 10-* - 10'' 


ocl (3) 


cxZ W 


(XZ (4) 


>107 


ocl (5) 


oc 1 (^) 


oc 1 (5) 



Note. — Given a set of ion fractions Xi, and a tem- 
perature T, the table shows how the cooling efficiency 
A{xi,T) depends on the gas metallicity Z. The entry 
"oc 1" means that the cooling efficiency is nearly inde- 
pendent of Z (but does depend on Xi and T). 

(1) FS line cooling. 

(2) H hja cooling. 

(3) H, He, and bremsstrahlung cooling. 

(4) Metal resonance line cooling. 

(5) Bremsstrahlung cooling. 

Lya emissions provide the gas cooling. For isochoric 
flows with Z > 1, flne-structure emissions dominate the 
cooling, and A is therefore proportional to 1/Z. For lower 
gas metallicities, and in _B = flows Lya emission dom- 
inates the cooling and A is independent of Z. 

Below ~ 7000 K, fine structure line emissions are the 
only available coolants we have included and A is there- 
fore proportional to Z for all gas metallicities. The de- 
pendence of A on Z is summarized in Table [6] below. 
These scalings will be useful in § [6l 

In Figure [201 we display the cooling efficiencies for Z 
between 10"'^ and 2, in the "hot radiative" and "non- 
equilibrium cooling" phases for the Tg = 5 x 10^ K (left 
panels) and Ts = 5 x 10^ K (right panels) shocks that 
we have considered, in the B — (lower panels) and 
strong-B limits (upper panels). We display results for 
temperatures between 5 x 10^ K, and the temperatures 
at which the photoabsorption plateau begins. 

Consider first the 5 x 10^ K shocks. For B = 
(isobaric) the gas compression reduces the photoioniza- 
tion parameter, and the non-equilibrium ion-fractions 
are mostly due to the recombination lags in the coUi- 
sionally ionized gas. In this limit, therefore, the cool- 
ing curves are almost identical to the non-equilibrium 
cooling curves presented in GS07 in which photoioniza- 
tion processes were excluded. However, as discussed in 
§131 in the strong-B (isochoric) limit, photoionization by 
the shock radiation significantly increases the ionization 
state of the gas above the non-equilibrium level that ob- 
tains in a purely radiatively cooling gas. More highly 
ionized species generally have more energetic resonance 
line transitions, and these are less efficiently excited by 
the thermal electrons at any given temperature (McCray 
1987). The cooling curves are therefore suppressed, com- 
pared to i? = shocks. 

For the hotter 5 x 10'' K shock, the shock radiation is 
more intense, and photoionization drives the ion fractions 
to still higher ionization states. Again, the more highly 
ionized gas cools less efficiently, and the cooling rates at 
a given temperature are therefore lower compared to the 
5 X 10^ K shocks. In the strong-B limit the intensity of 
the radiation is high enough that even minute amounts 
of neutral hydrogen allow for efficient photoabsorption 
heating, and the photoabsorption plateaus begin at a 
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t (s) T (K) 

Fig. 19. — The cooling (dark) and heating (gray) rates per unit volume (erg cm~^) for a strong-B 5 x 10® K at Z = 1. (a) Cooling 
and heating rates per unit volume as a function of time, assuming a post-shock hydrogen density no = 0.1 cm~^. (b) Cooling and heating 
rates per unit volume as a function of temperature. The dominant cooling and heating elements are indicated near the curves. 
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Fig. 20. — The cooling efficiencies, A (erg cm^ s as a function of temperature in the "hot radiative" and "non-equilibrium cooling" 
zones. Each panel shows results for gas metallicities between 10~^ and 2 times solar. For the range of temperatures displayed in this figure, 
heating is negligible. The curves end where the photoabsorption plateaus begin, (a) Tg = 5 x 10 K, for strong-B. (b) Tg = 5 x 10® K, for 
B = 0. (c) Ts = 5 X 10'^ K, for strong-B. (d) Ts = 5 x lO'' K, for B = 0. 
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higher temperature of 5 x 10 K. 

6. COOLING COLUMNS IN SHOCKED GAS 

In this section we present computations of the inte- 
grated metal-ion column densities that are produced in 
the post-shock cooling layers. In computing the cooling 
columns, we integrate from the shock front to a distance 
where the gas has cooled down to our termination tem- 
perature Tiow = 1000 K. 

As an example of what these integrations look like, in 
Figure [5T] we display the O^^ ion fraction and accumu- 
lated column density, in a strong-_B, 5 x 10^ K shock, 
for Z = 1. Panel (a) shows the O^"*" ion fraction as a 
function of time, for an assumed post-shock hydrogen 
density of 0.1 cm~^. Initially, the 0^+ abundance peaks 
as the gas is ionized after passing the shock front. This 
first peak is short, with a duration of less than 10^^ s. 
Later on, a second peak is formed as the gas cools be- 
hind the shock front. This second peak lasts for ~ 10^* s. 
Panel (b) shows the ion fraction as a function of tempera- 
ture. The first abundance peak occurs over such a narrow 
temperature range (T ~ Tg) that it cannot be resolved in 
panel (b). The second abundance peak is seen to occur 
between a temperature of ^ 10^ K, and ~ 10^ K. 

Panel (c) shows the accumulated 0^+ column density 
as a function of time. The initial abundance peak con- 
tributes an 0^+ column density of ~ 10^^ cm^^. The 
column density builds rapidly, and remains unchanged 
until > 10^^ s. The gray curve in panel (c) shows a scaled 
temperature as a function of time. During the hot radia- 
tive phase the accumulated column density remains un- 
altered, as higher ions dominate the oxygen abundance. 
The column remains constant until the non-equilibrium 
cooling phase, during which the second abundance peaks 
occurs. The column density then builds up rapidly to 
~ 10^^ cm~^. The 0^+ column keeps increasing in the 
photoabsorption plateau, until 0^+ finally recombines at 
T ~ lO** K. The accumulated column density as a func- 
tion of temperature is shown in panel (d). Panels (c) 
and (d) show that the contribution of the first abun- 
dance peak to the total column density in the post-shock 
cooling layers is negligible. 

More generally, the integrated metal ion column den- 
sity, Ni, through the post-shock cooling gas is given by. 



N,, 



with 



dt oc 



ni{Z Ac\Xivdt 



dT 

nA{xi,T, Z) 



(12) 



(13) 



In this expression, nn is the hydrogen density. Ad is the 
abundance of the element relative to hydrogen in a solar 
metallicity g the ion fraction, and v is the gas 

velocity. 

Since nnv is constant throughout the flow (eq. [3]), 
equations (fT2|) and (I13|l imply that 



dT 



(14) 



For strong-_B isochoric flows n — tiq everywhere, and 
Ni cx Z J XidT/A. However, for B — isobaric flows the 
density increases and the cooling times are shortened as 



the gas cools (see §[3]). Therefore, the column densities of 
ions produced in low temperature gas are substantially 
suppressed (by the shorter cooling times) in isobaric ver- 
sus isochoric cooling flows (e.g. Edgar & Chevalier 1986). 
For an ion that is produced at a characteristic tempera- 
ture T, the column density in a i? = shock is quenched 
by a factor approximately no/n{T) ~ T/Ts relative to a 
strong- i? shock. Because T/Tg is smaller for higher Tg, 
the quenching effect is more significant for higher shock 
velocities. 

Given an ion distribution Xi(T), equation (|14p implies 
that Ni is independent of if A cx Z (as occurs when 
metals dominate the cooling). When A (x Z, reduced 
elemental abundances are compensated by longer cool- 
ing times, and Z cancels out in equation (14). On the 
other hand, for ions that are abundant at temperatures 
where the cooling is dominated by hydrogen, helium, or 
bremsstrahlung emissions, equation l|14p implies that the 
column densities are proportional to Z. 

The actual behavior is more complicated because the 
ion distributions, Xi(T), are themselves a function of gas 
metallicity and shock temperature as discussed in § |4l 
For example, higher shock temperatures imply stronger 
photoionizing radiation fields, which affect the ionization 
states in the post-shock cooling layers. The gas metallic- 
ity affects the spectral energy distribution of the shock 
self-radiation, changing the ion fractions for different val- 
ues of Z . Departures from equilibrium ionization and the 
degree of over-ionization in the gas are also a function of 
gas metallicity. All these factors imply that the column 
densities depend on the metallicity even for temperatures 
where A cx Z. 

In Figures [22ll25l we display the full set of carbon-, ni- 
trogen, oxygen, and neon-ion column densities, for Z be- 
tween 10""^ and 2, for strong-B (upper panels) and B=0 
(lower panels) post-shock cooling layers. Left hand pan- 
els are for Tg = 5 x 10^ K, and the right hand panels 
for Ts = 5 X 10^ K. The shaded areas show ions that are 
still abundant at Tiow = 1000 K. The column densities 
of these ions (e.g. O and 0+) are sensitive to the choice 
of Tlow • The full set of column densities for all the metal 
ions that we consider are listed in Table [71 as described 
in Table H ^ 

Our integrated metal-ion cooling column densities are 
in good qualitative agreement with the recent results of 
Allen et al. (2008) for the case of a solar metallicity gas 
cooling behind a 600 km shock, in which B ~ 
(their model "J_nl_bO"). Remaining differences are likely 
due to differences in the assumed elemental abundances, 
and in the magnetic field strength. The computations 
presented in Allen et al. cover a different region of the 
parameter space, and do not show results for the other 
cases than we consider here (namely Vg = 2000 km s~^, 
sub-solar metallicities, and the limit of strong-i3). 

We now focus on the oxygen ion fractions shown in 
Figure [24] as an example. Table [6| predicts that the col- 
umn densities in low metallicity gas (< 10~^) are always 
proportional to Z except at very low temperatures. As 
can be seen from Figure [Ml the ratio between the oxygen 

The columns presented in Table [71 and Fsigures I22I25I assume 
that the shock is observed "face-on" (so that the line-of-sight is 
parallel to the shock velocity). These columns should be multi- 
plied by a geometrical correction factor, 1/cosd, for shocks that 
are viewed at an angle 9 to the normal direction. 
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t (s) T (K) 

Fig. 21. — O^^ ion fraction (upper panels) and accumulated column density (lower panels) as a function of time (left hand panels) and 
temperature (right hand panels), for a strong-B, T, = 5 X 10^ K, Z = 1 shock. The gray curve in panel (c) shows the temperature (rescaled 
by xlO^) as a function of time. 



TABLE 7 

Metal-ion Column Densities in a 5 x 10^ K, 2 times solar, 

STRONG-B-SHOCK 



ionization H He C N 

(cm~^) (cm~-^) (cm~^) (cm^^) 

I 5.7 X 10^0 4.4 X 10^'^ 8.8 x lO" 5.7 x 10^^ 

II 8.5 X lO^" 2.7 X 10^'-' 3.1 x 10^'^ 2.3 x 10^*^ 

III - 4.8 X 10^'-' 1.0 X 10^'^ 2.2 X lO^s 



Note. — The complete version of this table is in the electronic 
edition of the Journal. The printed edition contains only a sample. 
The full table lists metal-ion column densities for B = and in the 
strong-B limit, for shock temperatures of 5 X 10^ K and 5 X 10^ K, 
and for Z = 10~^, 10~^, 10~^, 1, and 2 times solar metallicity gas 
(for a guide, see Table |4j. 



ion column densities for Z = 10 ^ and the columns for 
Z = 10~^ is ~ 10 as expected. For example, in panel (a) 
the columns of 0^+, 0-*+, and 0^+ are - 8 x 10^-^, 
~ 7.5 X IQi'^, and - 7 X 10^'^ cm'^ for Z = 10"^ and 
~ 8 X 10^3^ ^ 7 5 X 1014^ and ~ 7 X 10" for Z = lO'^. 
The same applies for other shock temperatures and mag- 
netic field limits, as can be seen in panels (b)-(d). 

For Z > 1, Table [5] predicts that for a given ion dis- 
tribution Xi, the column densities are independent of Z. 
Figures [Mb andlMti show that for B = where photoion- 
ization is less significant, the column densities for Z — 1 
and for Z — 2 are indeed similar. Departures from equi- 
librium have only a very minor impact on the column 
densities. For example, in Figure [24h . the 0^^+ column 
is 9.6 X 10" cm-2 for Z = 2, and 9.1 x 10" cm'^ for 
Z = 1, due to the longer recombination lags for higher 
Z . However, in strong-_B isochoric shocks, photoioniza- 
tion introduces a Z dependence. For example, Figurc IMb 
shows that the 0^+ column density is 8.1 x 10^^ cm~^ 
for Z = 2, but is 3.1 x 10^^ cm-^ for Z = 1. Although 
A ^ Z, the 0^+ column density ratio of 2.6 is in this 
case even larger than the metallicity-ratio of 2.0. 



Figures [^^11^ clearly show the differences between 
strong-i? (isochoric) shocks and 5 = (approximately 
isobaric) shocks. For B — 0, the gas compression and 
shorter cooling times reduce the column densities of low 
and intermediate ions that are formed at T < Ts. A 
comparison of Figures [24k and [24h shows that while the 
column densities of the high ions (0®+, O^"*") is similar in 
both cases, the column densities of 0^+-0^+ are greatly 
suppressed in i? = models. For example, in solar metal- 
licity gas, the 0^+ column is, as expected, similar for the 
two cases. It is 5 x 10^^ cm~^ for a strong-B shock, and 
3 X 10^^ cm~^ for B — 0. However, for 0'^+ the column 
density is 8 x 10^^ cm~^ for strong-B, and 3 x 10^^ cm^-^ 
for B = 0, 3.5 orders of magnitude lower. The magnetic 
field therefore has a very strong impact on the column 
density ratios between high- and intermediate-ions, and 
X-ray and UV observations will be able to easily distin- 
guish between the cases of strong versus weak magnetic 
fields. 

The shock velocity affects the column densities pro- 
foundly. Equation[T2]implies that for ions that only exist 
below Ts, and given an ion distribution Xi{T), the inte- 
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Fig. 22. — Carbon ion-column densities in post-shock cooling layers, for Z from 10 ^ to 2 times solar, (a) Ts = 5 X 10® K, for strong-B. 
(b) Ts = 5 X 10'' K, for strong-S. (c) Ts = 5 X 10*^ K, for B = 0. (d) Ts = 5 X lO'^ K, for B = 0. The shaded areas show ions that are still 
abundant at Tio^ = 1000 K where we stop the computation. The column densities of these ions are sensitive to the choice of Ti^^. 



grated column density will be proportional to the shock 
velocity. However, since the shock self-radiation and the 
degree of photoionization in the gas depends on the shock 
temperature, this scaling is not exact. For i? = 0, pho- 
toionization is less important, and the ion distributions 
Xi{T) are therefore less sensitive to T,- For strong-S 
shocks, photoionization plays a major role in setting the 
ion fractions in the cooling gas, and the ion distributions 
are more affected by Tg. 

For example, panels (a) and (b) show that the column 
density of 0^+ is 8.1 x 10^^ cm ^ in a 5 X 10^ K shock of 
a two solar metallicity gas, compared to 2.2 x 10^^ cm~^ 
in a 5 X 10^ K shock. This is a factor of ^ 37 in the 
column density, compared to a factor of 10^/^ ^ 3 in the 
shock velocities. For = an additional factor T/Ts 
suppresses the ion column densities as discussed above, 
and the exact range of temperatures over which the ion 
is abundant will determine the column density ratio. 

7. METAL COLUMNS IN THE RADIATIVE 
PRECURSOR 

In this section we present equilibrium photoionization 
computations of the integrated metal-ion column densi- 
ties that are produced in the radiative precursors. The 
precursor gas is heated and photoionized by the shock 
self-radiation propagating upstream. The upstream ra- 
diation is gradually absorbed in the radiative precursor, 
producing a photoabsorption layer in which the heating 
and ionization rates gradually decline with (upstream) 



distance from the shock front. 

The ionization parameter in the radiative precursor is 
higher than in the post-shock photoabsorption plateau, 
because the gas density in the precursor is four times 
lower than the post-shock density (see 12. 2p . Typically, 
the temperature of the gas approaching the shock front 
is a few times 10'' K, but it may be as high as ^ lO'^ K 
for = 5 X 10^ K. In computing the precursor columns, 
we integrate from the shock front to a distance where the 
heating rate in the unshocked gas is sufficiently low that 
is reaches a temperature Tiow = 1000 K. 

We use Cloudy to construct equilibrium photoioniza- 
tion models for the radiative precursors. For the radia- 
tion flux entering the precursor, we use the post-shock 
flux at the position in which the post-shocked gas first 
becomes thick at the Lyman limit, ItMck, similarly to the 
scheme presented in § 12.51 For each shock model, we con- 
struct a precursor model assuming that a flux F^{lthick) 
enters a gas of density rto/4. The metal column densities 
are then Ni = J nuZAdXidl, where I is the upstream dis- 
tance from the shock front, and nn — no /A. In all cases, 
the precursors are assumed to have a constant density. 

Figures [26ll29l display the full set of carbon-, nitrogen-, 
oxygen-, and neon- ion column densities, for Z between 
10""^ and 2, for the radiative precursors of strong- i? (up- 
per panels) and B — (lower panels) shocks. Left hand 
panels are for Tg ~ 5 x 10^ K, and right hand panels are 
for Ts — 5 X 10^ K. The shaded areas show ions that are 
still abundant at Tiow = 1000 K. The column densities 
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of these ions (e.g. O and 0+) are sensitive to the choice 
of Tiow The full set of column densities for all the metal 
ions that we consider are listed in Table [51 as described 
in Table [4] Our integrated precursor column densities 
are in good qualitative agreement with the recent results 
of Allen et al. (2008) for the case of a solar metallicity 
gas cooling behind a 600 km s^^ shock, in which B ~ 
(their model "J_nl_bO"). 

Figures [26ll29l show that significant column densities 
are created in the radiative precursors. For the high- 
ions, which coUisional abundance peak at T ~ Tj, the 
column densities are dominated by the post-shock cool- 
ing layers. However, the higher ionization parameters 
in the lower density precursors allow for a very efhcient 
production of lower ions that are created below Tg. The 
precursor columns often exceed the post-shock cooling 
columns, especially for ions that are maintained by pho- 
toionization in the post-shock gas. For example, in solar 
metallicity gas, the ratio of post-shock to precursor col- 
umn density of the high ion O VIII, is iVpostZ-^i 



* pre 



1000 

for a = 5 X 10^ K strong-B shock. The O VIII column 
is completely dominated by the post-shock gas. How- 
ever, for O IV this ratio is 0.25, so that most of the O IV 
column is created in the precursor. 

Higher shocks velocities produce harder and more in- 
tense radiation-fields, resulting in efficient production of 
high-ions and in deeper, higher column-densities, radia- 
tive precursors. For example the upper panels in Fig- 
ure [28] show that when Z = 1, the O VIII column is 
8.5 X lO^^ cm-2 for T, = 5 x lO"^ K, but 7.7 x lO^^ ^^-2 



for Ts^5x 10^ K. The O IV column is 5.9 x lO" cm^^ 
forT^ = 5x 10^ K, but 3.1x101^ cm^^ ^^^j,^ ^ 5x106K. 

The column densities produced in the radiative precur- 
sors of strong- i? shocks are similar to those produced in 
B — precursor. The radiative fluxes in the two cases 
differ only by a factor 5/3 (see Figure [7]). A compari- 
son of the upper and lower panels in Figures [^M^ shows 
that the precursor columns are indeed similar. 

For a given radiation field and absorption depth, the 
metal-ion column densities scale linearly with gas metal- 
licity. However, since both the spectral energy distribu- 
tion and the gas opacities are functions of the gas metal- 
licity, some deviations from this linear correlation are 
expected. For example, lower-Z shocks produce harder 
radiation fields, which are more efficient in producing 
high ions. Figures [26ll29l shows that the metal-ion col- 
umn densities generally scale with gas metallicity, but 
departures from exact linear correlation are apparent in 
some cases. 

While the column densities in the radiative precursor 
may far exceed the columns produced in the post-shock 
cooling layers, kinematic differences may distinguish the 
two components. For high-ions created at Tg, the 
kinematic offset between the radiative precursor and the 
post-shocked gas is 3us/4, which is of order the ther- 
mal width, fons produced at lower temperatures have 
narrower widths, making it possible to distinguish the 
precursor and post-shock components. In i? = shocks, 
the post-shock low-ion columns are greatly suppressed 
(see § [S]), and are orders of magnitudes lower than the 
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TABLE 8 

Precursor Column Densities in a 5 x 10^ K, 2 times solar, 

STRONG-B-SHOCK 



ionization 


H 


He 


C 


N 




(cm-2) 


(cm~^) 


(cm~^) 


{cm~^) 


I 


2.2 X 10^1 


1.7 X 10^0 


7.8 X 10" 


2.4 X 10^'' ... 


II 


2.3 X 10^1 


1.5 X 10^0 


1.2 X 10^** 


8.4 X 10^^ ... 


III 




5.5 X 10^^ 


5.2 X 10^^ 


9.1 X 10^'^ ... 



Note. — The complete version of this table is in the electronic 
edition of the Journal. The printed edition contains only a sample. 
The full table lists precursorn column densities for B = and in the 
strong-B limit, for shock temperatures of 5 X 10^ K and 5 X 10^ K, 
and for Z = 10~^, 10~^, 10~^, 1, and 2 times solar metallicity gas 
(for a guide, see Table |4)l. 

precursor columns. In strong-i? isochoric shocks, signif- 
icant post-shock low-ion columns are expected to form 
with a velocity offset of 3ws/4 relative to the precursor, 
and with the same velocity centroid as the post-shock 
high-ions. 

8. SUMMARY 

In this paper we present new computations of the 
metal-ion column-densities produced in post-shock cool- 
ing layers behind fast, radiative shocks. We have con- 
structed a new (one-dimensional, steady) shock code 
in which we explicitly follow the non-equilibrium ion- 
ization states of post-shock gas containing H, He, C, 
N, O, Ne, Mg, Si, S, and Fe. We present results 



for initial post-shock temperatures, Tg, of 5 x 10® and 
5 X 10^ K, corresponding to shock velocities, Vs, of 600 
and ~ 2000 km s^^. We consider shocks in which there is 
no magnetic field {B — 0) for which the cooling flows are 
approximately isobaric (Poo = (4/3)Po)- We also con- 
sider shocks in which the magnetic field dominates the 
pressure everywhere ("strong-P" ) for which the flows are 
isochoric (rtoo = "o)- We assume that the gas is dust free, 
and we present results for metallicities Z ranging from 
10~^ to twice the solar photospheric abundances of the 
heavy elements. 

For the shock temperatures that we consider, the post- 
shock gas emits energetic radiation that may later be 
absorbed by cooler gas further downstream, and by the 
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Fig. 25. — Same as Figure [22] but for Neon. 



unperturbed gas that is approaching the shock front. 
The shock self-radiation significantly affects the ioniza- 
tion states and thermal properties of the gas. Its absorp- 
tion by the downstream gas creates a photoabsorption 
zone, in which the temperature and ionization states of 
the gas are set by the shock self-radiation. 

In following the time-dependent ion fractions, we rely 
on the code and work presented in Gnat & Sternberg 
(2007). In this paper, in addition to calculating the non- 
equilibrium ionization and cooling, we also follow the ra- 
diative transfer of the shock self-radiation through the 
post-shock cooling layers; take into account the resulting 
photoionization and heating rates; follow the dynamics 
of the cooling gas; and self-consistently compute the ini- 
tial photoionization states in the precursor gas. We use 
up-to-date rate coefficients for all of the atomic ioniza- 
tion and recombination processes, and for the energy loss 
and absorption mechanisms. The equations and our nu- 
merical method are presented in § [2l 

In § |3] we discuss the shock structure and emitted ra- 
diation, and we discuss how these quantities depend on 
the controlling parameters, including the gas metallic- 
ity, shock velocity, magnetic field and gas density (see, 
e.g. Fig. 1). For B ^ shocks (nearly isobaric) the gas 
is compressed and decelerated as it cools. Because the 
cooling time-scale is proportional to 1/n, the final evolu- 
tionary stages of isobaric shocks occur much more rapidly 
than in " strong- i?" isochoric shocks. The effects of pho- 
toionization by the shock radiation are much more impor- 
tant in isochoric shocks because the ionization parameter 



remains large in the downstream absorbing layers. 

The evolution of the post-shock cooling layers is sig- 
nificantly affected by the gas metallicity in several ways. 
High metallicity shocks cool more rapidly, because for 
10* ^ T < 10^ K the cooling efficiency is dominated 
by metal-line cooling. The degree to which recombina- 
tion lags behind cooling, and the resulting departures 
from ionization equilibrium therefore depend on the gas 
metallicity. Second, while the total energy emitted by the 
cooling gas is independent of Z, and equals the energy 
flux entering the shock, the spectral energy distribution 
of the shock self-radiation is a function of gas metallicity. 
For high-Z, a high fraction of the input energy is radiated 
as line-emission, while for \ow-Z, the relative contribu- 
tion of lines is small, and most of the initial energy flux 
is radiated as thermal bremsstrahlung continuum. At 
low Z the shock self-radiation is harder. The photoion- 
ization states of the downstream layers depend on the 
metallicity-dependent spectral energy distributions. 

In § [¥] we present our computations of the non- 
equilibrium ionization states in the post-shock cooling 
layers. In Tabic [3] wc list the time-dependent ion frac- 
tions as a function of time and temperature for the vari- 
ous shock models that we consider. As the gas enters the 
shock, its ionization state rapidly adjusts toward CIE at 
the shock temperature. Later the gas gradually cools, re- 
combines, and radiates its thermal energy. We describe 
how photoionization by the shock self-radiation affects 
the ion fractions in the cooling gas. We find that in 
strong- i? (isochoric) shocks, photoionization plays a cen- 
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Fig. 26. — Carbon ion-column densities in radiative precursor, for Z from 10 to 2 times solar, (a) Ts = 5 X 10^ K, for strong-B. 
(b) = 5 X 10^ K, for strong-B. (c) Ts = 5 X 10*^ K, for B = 0. (d) Ts = 5 X lO'^ K, for 5 = 0. The shaded areas show ions that are still 
abundant at T\q.^ = 1000 K where we stop the integration. The column densities of these ions are sensitive to the choice of Tio„. 



tral role in setting the ionization states in the coohng 
layers, even when the gas is still optically thin. Pho- 
toionization increases the abundances of high-ions at low 
temperatures, above the enhancements due to the re- 
combination lags. In those phases where photoionization 
dominates, the ionization states remain close to pho- 
toionization equilibrium. For B — Q (isobaric) shocks, 
the gas compression significantly suppresses the role of 
photoionization. 

We further demonstrate how departures from ioniza- 
tion equilibrium alter the gas ionization. We find that for 
strong-i? shocks, departures from photoionization equi- 
librium occur at temperatures between 2 x 10'* K and 
^ 10^ K, but the ion fractions differ by < 25%. However, 
when 5 = 0, departures from equilibrium ionization are 
larger due to the smaller effect of photoionization, and 
the non-equilibrium ion-fractions are mostly due to the 
recombination lags in coUisionally ionized gas. We also 
study the how the ion distributions depend on the gas 
metallicity, shock velocity, and magnetic field. 

We present our computations of the radiative cooling 
and heating efficiencies in § [5l We list our results in 
Table [5] and display them in Figure [20l We discuss the 
different cooling and heating processes that operate in 
the various cooling "zones" (or "phases") and shown in 
Figures [T] and ID 

Finally, in § [6] and § [7] we present our results for the 
integrated metal- ion column-densities. In §[6] we discuss 
the columns that are produced in the post-shock cool- 



ing gas. We give results for gas metallicities between 
10"'^ and 2 times the solar abundances. The computed 
column densities are listed in Table [T] We show that 
the cooling column densities created in the post shock 
cooling layers are strong functions of the magnetic field 
intensity, shock velocity and gas metallicity. Ionic ra- 
tios are useful as diagnostic probes. The predicted col- 
umn densities show that UV and X-ray absorption-line 
observations of shocked gas may be used to probe the 
intensity of magnetic fields in interstellar and intergalac- 
tic shocks, and to infer the physical properties of gas in 
fast radiative shock waves. In § [7] we compute the inte- 
grated equilibrium metal-ion column densities produced 
in the upstream radiative precursors. Our precursor col- 
umn densities are listed in Table [5] For "low-ions" , the 
column densities produced in the radiative precursors are 
often comparable to or greater than the post-shock cool- 
ing columns. However, kinematic differences between 
the precursor and post-shock components may provide 
a means to distinguish observationally between the two 
contributions. 
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Fig. 27. — Same as Figure [22] but for Nitrogen. 
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Fig. 29. — Same as Figure [22] but for Neon. 



